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FOREWORD 


This  document  is  Volume  III  of  a  four  volume  report. 
The  four  volume  report  is  titled  MNAVSTAR  Global 
Positioning  Systems  Special  Studies  and  Engineering 
Program",  Volume  III  presents  the  results  of  two 
studies:  One  titled,  "Direct  Acquisition  Study" 
previously  submitted  on  11  July  1975  at  STI-TR  7115-2 
the  other  titled  "Manpack  Study"  previously  submitted 
on  11  July  1975  as  STI-TR  7115-3. 
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GPS  SPECIAL  STUDIES  AND  ENGINEERING  PROGRAM 


TASK  II  -  DIRECT  ACQUISITION  STUDY 


1.0  INTRODUCTION 

The  purpose  of  the  direct  acquisition  is  to  determine  the 
feasibility  of  a  user  acquiring  a  P  signal  directly  (i.e, 
without  the  aid  of  the  C/A  signal)  in  the  presence  of 
jamming. 

This  investigation  was  approached  in  the  following  manner: 

(1)  Determine  direct  acquisition  performance 
based  on  presently  anticipated  system  elements 
including  the  constellation  selection  program, 
satellite  oscillators  as  specified,  reasonable 
user  oscillators,  etc. 

(2)  Conduct  parametric  investigations  to  determine 
requirements  on  critical  system  elements  for 
acceptable  direct  acquisition  performance. 

This  work  was  completed  during  August,  1974  and  reported  in 
STI  memos  GPS-008,  6PS-014,  and  GPS-018,  from  which  this 
report  is  derived. 
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Conclusions  and  Recommendations 


•  The  receiver  local  oscillator  is  the  limiting 

factor  in  performing  direct  acquisition.  A 

receiver  local  oscillator  with  overall  stability 
-9 

of  10  or  better  appears  to  be  required. 

•  The  capability  to  remotely  update  the  user  clock 
in  phase  and  frequency  is  required.  This  appears 
feasible. 

•  Direct  acquisition  appears  feasible  for  the  user 
with  velocity  uncertainty  <5  MPH  and  with  an 
oscillator  which  meets  the  requirements  stated 
above.  The  low  velocity  uncertainty  user 
employing  a  parallel  receiver  with  oscillator 

_9 

stability  <10  is  expected  to  acquire  all  four 
signals  within  about  300  seconds  with  probability 
^.9  twenty-four  hours  after  clock  synchronization. 
Direct  acquisition  is  considerably  more  difficult 
for  the  user  with  velocity  uncertainty  as  great 
as  60  MPH.  In  this  case,  an  oscillator  with 
stability  on  the  order  of  10~^  is  required  if 
direct  acquisition  is  to  be  practical  twenty- 
four  hours  after  clock  synchronization. 

Direct  acquisition  within  25  hours  of  an  accurate 
clock  update  appears  to  be  feasible  for  the  manpack 
receiver.  However,  the  oscillator  requirements  and 
addition  receiver  complexity  composed  are  such  that 
this  mode  of  operation  is  not  recommended  for  the 
manpack  user  equipment  to  be  developed  in  the  near 
future . 
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2.0  User  Baseline  and  Parameters 

Since  the  acquisition  time  results  are  dependent  on  both 
the  assumed  system  and  scenario,  it  is  important  that  baselines 
be  defined.  It  seems  that  direct  acquisition  would  be  desirable 
in  two  representative  situations: 

(1)  The  user  is  not  navigating  and  wishes  to  enter 
the  system  in  the  presence  of  jamming 

(2)  The  user  is  navigating  but  wishes  to  acquire 
a  P  signal  in  the  presence  of  jamming  due  to 
satellite  handover,  monentary  lost-of-lock,  etc. 

Obviously  the  user  in  the  second  situation  has  much  more  current 
and  accurate  system  information  than  the  user  in  the  first 
situation  and  consequently  a  much  less  difficult  reacquisition 
problem.  The  second  situation  will,  for  the  moment,  be 
considered  as  a  special  case  of  the  first  situation  and  not 
be  discussed  further  at  this  time.  We  will  treat  the  case 
where  the  user  is  not  navigating  but  wishes  to  enter  the  system 
in  the  following  material. 

2.1  Direct  Acquisition  Procedure 

The  user  does  not  know  the  exact  code  phase  of  the  received 
signal  and  must  therefore  perform  a  code  search  by  stepping 
the  phase  of  the  reference  code  with  respect  to  the  incoming 
code  until  correlation  is  achieved.  This  search  process 
is  relatively  slow  in  the  presence  of  jamming  (<10  chips/sec 
for  significant  jamming  levels),  requiring  that  the  initial 
code  phase  uncertainty  be  held  to  a  minimum  for  reasonable 
acquisition  times.  For  example,  a  3000  chip  uncertainty, 
corresponding  to  an  uncertainty  of  +150  fi8  in  code  timing 
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searched  at  10  chips /sec,  would  yield  a  maximum  acquisition 
time  of  300  sec.  This  time  uncertainty  probably  represents 
an  upper  bound  on  what  is  acceptable.  The  user  experiences 
uncertainty  in  code  arrival  time  due  to  uncertainties  in: 

•  system  time 

•  satellite  time 

•  satellite/user  range 

•  ionospheric  delay 

It  is  obvious  from  the  above  discussion  that  the  user  must 
have  a  fairly  accurate  estimate  of  system  time  (probably 
with  an  error  somewhat  less  than  150  /xs).  Thus  we  may  assume 
that  the  user  will  require  a  synchronized  clock  as  well  as 
some  means  for  estimating  satellite/user  range. 

Further,  the  code  search  may  be  carried  out  with  only  a 
relatively  small  uncertainty  in  carrier  frequency  (around 
100  or  200  Hz).  Larger  uncertainties  require  parallel 
correlators  or  additional  search  time.  Thus  the  user  is 
required  to  estimate  satellite/user  doppler.  The  user  also 
requires  an  accurate  estimate  of  received  code  frequency 
(within  1  or  2  Hz)  in  order  to  achieve  the  desired  code 
search  rate.  This  estimate  is  derived  along  with  the  carrier 
f r'  *iency  estimate. 

The  discussion  that  follows  will  assume  the  baseline: 

•  The  user  must  incorporate  a  clock  which  requires 
synchronization 

•  Some  time  elapses  between  clock  synchronization 
and  the  direct  signal  acquisition  attempt 
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The  user  employs  his  constellation  selection 
program  to  estimate  satellite /user  range  and 
range-rate 

The  satellite  clock  performance  is  consistent 
with  the  present  satellite  specification 

The  user  employs  a  single  noncoherent  correlator 
with  input  bandwidth  Bj. 


3.0  CLOCK  SYNCHRONIZATION 

Two  types  of  user  clock  synchronization  may  be  considered.  The 
first  is  "direct"  synchronization  from  a  local  standard , 
such  as  a  cessium  beam  clock  located  at  an  equipment  depot, 
airport,  on  board  ship,  etc.  The  synchronization,  in  this  case, 
is  readily  accomplished  to  within  the  accruacy  of  the  local 
standard  and  requires  very  little  complexity  on  the  part  of  the 
user. 

Unfortunately,  direct  synchronization  at  relatively  frequent 
intervals  may  be  inconvenient  or  impossible  for  many  users. 

The  user,  when  navigating,  may  perform  "remote"  synchronization 
of  his  clock  as  follows: 

•  The  user  receives  satellite  time  from  the 
P  signal  data  corresponding  to  code  time 
ticks  every  1.5  sec. 

•  Since  the  user  knows  his  own  position  and  the 
satellite  orbit,  he  may  compute  propagation 
delay  time. 

•  The  user  also  receives  satellite  clock  phase 
correction  coefficients  to  allow  calibration 
with  system  time. 

Thus  the  navigating  user  may  set  his  clock  to  within  a  few 
qs  of  correct  system  time  using  the  signal  from  any  satellite. 
Further,  he  may  make  a  frequency  correction  by  synchronizing 
with  the  received  carrier,  correcting  for  user /satellite  doppler, 
and  correcting  for  satellite  clock  frequency  error  (received  on 
the  P  signal  data)  •  A  block  diagram  of  the  clock  synchroniza¬ 
tion  is  shown  in  Figure  1. 


Figure  1  Clock  Synchronization 


.... _ 


The  above  synchronization  procedure  can  actually  be  accom¬ 
plished  when  only  one  satellite  signal  is  being  received, 
however  the  synchronization  accuracy  is  degraded  due  to  the 
uncertainty  in  user  position  and  velocity. 

In  summary,  remote  snychronization  when  the  user  is  navigating 
can  be  very  accurate,  with  negligible  error  in  phase,  and 
frequency  error  limited  primarily  by  the  accuracy  with  which 
the  user  oscillator  may  be  set. 
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SIGNAL  ACQUISITION 


Let  us  assume  that  after  clock  synchronization,  the  user 
receiver  is  turned  off,  except  the  clock,  for  some  period 

of  time  and  the  user  then  attempts  direct  P  signal  acquisi¬ 

tion.  The  user  must  first  estimate  received  code  phase  based 
on  estimates  of 

•  system  time 

•  satellite/user  range. 

The  code  phase  estimate  will  be  in  error  due  to: 

•  error  in  satellite  position  estimate 

•  error  in  user  position  estimate 

•  error  in  user  estimate  of  system  time 

•  error  in  satellite  clock  with  respect 

to  system  time 

•  ionospheric  delay 

These  error  sources  are  discussed  below. 

4.1  Satellite  Position  Estimate 

Let  us  assure  that  the  user  relies  on  his  initial  orbit 
determination  model  for  a  satellite  position  estimate.  The 
accuracy  of  such  a  model,  employing  the  Kepler  two-body 
equations,  is  discussed  in  Reference  (1).  The  use  of  average 
values  for  orbit  elements  appears  to  be  the  best  approach, 
since  it  gives  a  relatively  stable  value  of  maximum  satellite 
position  error  for  12  to  16  days  after  element  up-date.  The 
in-track  position  error  varies  from  about  40,000  feet  to 
70,000  feet  over  a  sixteen  day  period.  The  change  in  satellite/ 
user  range  due  to  in- track  satellite  position  errors  may  be 
computed  by  employing  the  equations  (Ref.  (1)). 


R.E.  Orrf  TRW,  DNSDP  -  RED-155 
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2  2  2 

r  83  r  +  r  -  2r  r  cos  a 
e  s  e  s 


and  A  ©  «  ^2. 

rs 


where 


rg  is  the  satellite  geocentric  radius,  87.2  x  10  6  ft 
re  is  the  earth’s  radius,  20.9  x  10^  ft 

r  is  the  range  from  user  to  satellite 

AP  is  the  variation  in  satellite  position  along 
the  direction  of  motion. 


Satellite 


E  Earth 
Center  ° 


User 


Differentiation  Eq  (1)  with  respect  to  a  gives 


dr  _  rers 
cla  r 


sin  « 


and  substituting  Eq  (2)  yields 
UP 

AR  ■  -  sin  o 

Substituting  numerical  values. 


AR  as  .23AP  when  <»  “  r/2 
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Thus  we  anticipate  errors  in  the  range  estimate  on  the  order 
of  ±14,000  ft,  based  on  an  error  of  ±60,000  ft  in  the 
satellite  position  computation  (i.e.  a  time  uncertainty  of 
±14  fi s). 
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System  Time  Estimate 


The  user  clock  will  be  in  error  due  to 

•  initial  time  synchronization  errors 

•  initial  frequency  synchronization  errors 

•  user  clock  drift. 

The  initial  time  synchronization  errors  are  expected  to  be 
negligible  as  discussed  in  Sec.  3. 

An  initial  frequency  synchronziation  error  may  exist  due  to 
errors  in  the  doppler  and  satellite  oscillator  correction, 
as  well  as  an  error  in  setting  the  oscillator.  The  doppler 
estimate  and  satellite  frequency  corrections  should  introduce 
a  maximum  combined  error  of  about  ±2  x  10~^,  based  on  the 
requirement  for  range-rate  accuracy  of  .06M/sec  for  the  Class 
X  and  Y  user  (the  Class  Z  user  accuracy  is  not  yet  specified). 
An  oscillator  setting  accuracy  of  ±10"^  will  be  assumed  (See 
Appendix  A). 


Once  the  user  oscillator  has  been  set,  there  will  be  frequency 
variations  (for  quartz  crystal  oscillators)  generally  categorized 
as  being  due  to  temperature  variation,  voltage  variation,  load 
variation,  short  term  instability,  and  aging. 
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For  the  purpose  of  our  investigation,  let  us  assume  that  the 
user  oscillator  frequency  is  given  by 


fQ  +  Af  +  at 


where 


fQ  =  nominal  system  standard  frequency 

Af  *  a  constant  user  frequency  bias 
a  *  frequency  aging  rate. 


The  frequency  bias,  Af,  will  be  considered  to  be  a  Gaussian 
random  variable  with  standard  deviation 


‘r2  +  ‘v2  +  ‘i  +  ‘sis  +  V 


where 


temperature  stability 
voltage  stability 


load  stability 


short  term  stability 


setability 


A  discussion  of  parameters  assumed  for  the  above  error  sources  is 
presented  in  Appendix  A*  Three  baseline  reference  oscillators 
are  considered  for  this  investigation: 


A  low  power  crystal  oscillator  (LPO)  with  fast 
warmup  and  proportional  oven  temperature  control 
for  which  <jKC  -o 


6  x  10*®  and  jr- 


+  10"®/day. 
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2 .  An  ultra  stable  crystal  oscillator  (USO) ,  with 
a  double  proportional  control  oven,  for  which 

jr^  ■  1.4  x  10  ^  and  ^  m  +  10  ^/day. 

£o  ro 

.n 

3.  An  atomic  standard  (AS)  with  accuracy  of  10  . 

These  baseline  oscillators  are  discussed  in  more  detail  in 
Appendix  A. 

4.3  User  Position  Estimate 

The  error  in  the  user  position  estimate  will  depend  on  individual 
circumstances.  Position  uncertainty,  for  present  purposes,  will 
be  considered  to  be  a  circle  with  5  n.mi  radius,  consistent  with 
Table  II  of  the  User  System  Segment  Specification,  SS-US-101A. 

This  will  yield  a  worst  case  timing  uncertainty  of  approximately 
+  30  /is. 

4.4  Satellite  Clock  Error 

The  error  in  the  satellite  clock  will  be  less  than  +  10  /is, 
according  to  present  system  specifications. 


[onosi 


The  ionospheric  delay  is  a  few  chips  at  most  and  will  therefore 
be  neglected. 


The  anticipated  time  uncertainties  are  summarised  in  Table  1. 
There  is  some  question  as  to  how  these  time  uncertainties  should 
be  combined  to  give  the  most  meaningful  results.  For  present 
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purposes,  let  us  assume  that  each  error  is  a  Gaussian  random 
variable  with  2*  value  equal  to  the  maximum  value  shown  in  Table  1 
(with  the  exception  of  the  user  clock  error) •  Hie  total  time 
uncertainty  will  also  be  considered  to  be  Gaussian  with  standard 
deviation 

at  •  JWFm  ♦  y  * 

where  (UC)  is  the  variance  of  the  appropriate  user  clock. 

Using  the  numerical  values  of  Table  lt  we  have 

„aT-  \J 299  +  'T2  (UC) 

This  quantity  is  plotted  vs  time  in  Figure  2  for  the  different 
user  oscillators. 


j&aaa  fcar- -seat . <•,  ^  — 


4-7 

■a“atol — — 


— 


4 . 7  Carrier  Frequency  Estimate 


The  rate  the  user  can  search  the  received  code  timing  uncertainty, 
depends  partly  on  the  accuracy  with  which  the  signal  carrier 
frequency  is  known.  The  received  carrier  frequency  may  be  off 
nominal  due  to  satellite  oscillator  frequency  offset  and  satellite/ 
user  doppler.  Further,  frequency  error  in  the  user  oscillator 
will  also  appear  as  a  signal  offset.  The  problem  of  estimating 
signal  carrier  frequency  is  essentially  the  same  as  that  encountered 
for  C/A  signal  acquisition  and  is  discussed  in  Reference  (2) . 

These  estimates,  presented  in  Reference  (2),  may  be  used  when 
considering  direct  signal  acquisition,  except  that  the  error  due 
to  the  user  oscillator  is  anticipated  to  be  reduced  by  the 
synchronization  process  as  described  previously. 


A  summary  of  estimated  frequency  errors  for  the  direct  acquisi¬ 
tion  problem  is  shown  in  Table  2.  These  error  sources  will  be 
handled  in  the  same  fashion  as  the  time  error  sources,  i.e.  the 
errors  will  be  assumed  Gaussian  with  a  2*  value  equal  to  the 
maximum  value  shown.  The  standard  deviation  of  the  frequency 
error  will  be  expressed  as 


f^Y 


Using  the  values  in  Table  2  gives 


m 


rf  (UC) 


6364  +  (UC)  dynamic  user 


464  +  <y(iJC) 


stationary  user 


where  Of(UC)  is  the  one  sigma  frequency  error  of  the  user  clock. 
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4.8  Code  Frequency  Estimate 

The  rate  at  which  the  user  searches  the  received  code  may  be 
controlled  only  to  within  the  accuracy  with  which  the  received 
and  nominal  code  rates  are  known,  i.e.  the  user  must  reduce 
the  nominal  search  rate  below  the  desired  maximum  search  rate 
by  the  amount  of  uncertainty  in  the  code  frequency  estimate. 
Hie  factors  involved  in  estimating  code  frequency  are  the  same 
as  those  discussed  in  connection  with  estimating  carrier 
frequency.  Hie  accuracy  of  the  code  frequency  estimate  may  be 
determined  by  multiplying  the  carrier  frequency  uncertainty  by 
the  ratio  of  the  code  to  carrier  frequencies,  ■j-jg'  .  Thus  we 
expect  a  maximum  error  in  the  code  frequency  estimate  of 
about  +  .5  Hz. 
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Clock  Svchronization  Summary 


The  user  must  have  some  means  for  synchronizing 
an  internal  clock.  It  is  feasible  to  perfrom 
clock  synchronization  with  negligible  error 
during  navigation. 

The  user  will  have  a  1  or  timing  error  ~17  i*s 
(170  chips)  in  estimating  receive  code  timing, 
assuming  no  user  clock  error. 


User  induced  doppler  is  a  major  source  of 
frequency  uncertainty. 


5.0  ACQUISITION  TIMES  FOR  DIRECT  P-SEQUENCE  ACQUISITION 

Direct  acquisition  of  the  F-sequency  by  a  non -navigating  GPS 
user  is  discussed  in  some  detail  in  the  preceeding  sections. 

The  rate  at  which  the  received  code,  or  time  uncertainty, 
may  be  searched  is  derived  in  App.  B  with  probability  of  detec¬ 
tion  (PD)  IF  bandwidth  (Bj),  SNRj,  and  probability  of  false 
alarm  (Pp^)  as  parameters.  The  correlation  detector  in  Figure 
3  is  assumed. 


Figure  3  In-Lock  Detector 
5.1  Acquisition  Time 

The  standard  deviation  of  the  received  signal  code  time 
uncertainty  is  shown  in  Figure  2  vs  time  from  clock  synch¬ 
ronization  for  various  user  oscillator  stabilities.  Assuming 
that  the  time  error  is  a  Guassian  random  variable,  the  search 
interval  of  ±2<rAT  will  include  the  correct  time  position  with 
probability  .977. 

The  rate  at  whicl\  the  code  may  be  searched  is  plotted  vs  C/Nq 
with  Bj  as  a  parameter,  in  Figure  4  for  P^  ■  .94  and  Pp^  *  .005. 
A  correlator  loss  of  1  dB  is  assumed.  The  required  IF  bandwidth 
depends,  of  course,  on  the  carrier  frequency  uncertainty  (given 
by  Eq.  2)  and  is  taken  to  be  Bj  %  4 in  the  results  that  follow. 


The  time  to  acquire  one  signal  may  now  be  determined  for 
a  given  C/Nq  (or  J/S)  by  finding  the  corresponding  search 
rate  from  Figure  3  where  Bj  is  selected  based  on  (Eq.  2), 
and  the  time  interval  to  be  searched  is  given  by  Eq.  1. 

In  the  results  that  follow,  we  choose  a  time  search  interval 
of  +  2c an  IF  bandwidth  of  +  2 and  a  s  .94.  Thus 
the  probability  of  acquiring  the  signal  in  a  single  pass  is 
0.9.  The  results  plotted  are  maximum  acquisition  time,  i.e. 
time  to  search  the  interval  -  2c^T  to  +  2o^T  (the  increase  in 
acquisition  time  due  to  false  alarms  is  assumed  negligible) . 

The  optimum  search  pattern  would  start  at  the  middle  of  the 
time  uncertainty  interval  and  spiral  outward.  However,  it  is 
considerably  more  straightforward  to  implement  a  search  that 
starts  at  one  end  of  the  uncertainty  interval  and  searches  in 
one  direction.  In  this  case,  the  average  acquisition  is  one- 
half  of  the  maximum  time. 

Maximum  acquisition  time,  Tacq(max),  is  shown  vs  time  from 

clock  synchronization  in  Figure  5  for  a  user  with  velocity 

uncertainty  <  5  MPH  (corresponding  to  a  large  segment  of  the 

manpack  population) .  Oscillator  stability  is  seen  to  be  the 

overriding  factor  in  determining  performance,  with  a  stability 

.0 

of  10  or  better  being  required  for  acceptable  performance. 

Similar  curves  are  shown  in  Figure  6  for  the  user  with  a 
maximum  uncorrected  doppler  of  +  160  Hz  (corresponding  to  a 
velocity  uncertainty  of  30  m/sec  horizontal  and  6  m/sec 
vertical  or  60  MPH.  The  search  rate  of  5  chips/sec  is 
determined  as  follows.  An  IF  bandwidth  -  350  Hz  is  required  to 
accommodate  the  total  frequency  uncertainty.  Interpolating 
between  the  curves  of  Figure  4  we  find  a  maximum  search  rate 
of  7  chips /sec  at  j/S  -  45  dB  (C/Nq  -  27  dB-Hz) .  However,  a 
code  frequency  uncertainty  -  ±  1  chip /sec  exists  due  to  the  same 
sources  that  cause  an  uncertainty  in  the  carrier  frequency. 
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Figure  5  Acquisition  Time  for  User  with 

High  Velocity  Uncertainty  ( <  60  MPH) 


J/S 

=  45  dB 

(C/N0 

=  27  dB-Hz 

BI 

=  350  Hz 

“  6  Chips /sec 

PD 

■  .94,  PpA  -  .0 

1  dB  Correlation  Loss 

.cl 


*  *****  *f**35^S ^ 


Thus  the  nominal  search  rate  is  reduced  to  6  chips/sec  to 
ensure  that  the  maximum  search  rate  is  not  exceeded.  The 
results  indicate  that  relatively  large  uncorrected  doppler 
causes  a  drastic  increase  in  acquisition  time  (by  a  factor 
~2.8)  and  limits  useful  performance  (Tacq  -  3  min)  to  about 
five  hours  from  clock  synchronization  for  an  oscillator  with 
10~9  stability. 
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Acquisition  of  Additional  Signals 


Certain  user  receivers  may  attempt  to  acquire  four  satellites 
in  parallel,  in  which  case  the  maximum  time  to  acquire  four 
satellites  is  the  same  as  the  maximum  time  to  acquire  one 
satellite  (although  the  probability  of  acquiring  all  four  on 
the  first  attempt  is  only  about  0.65). 

In  case  the  receiver  does  not  acquire  all  four  signals  (or 
if  the  user  employs  a  sequential  receiver)  the  acquisition  of 
a  single  satellite  signal  allows  the  user  to  update  his  clock 
as  discussed  in  Section  3.  The  other  user  may  still 
experience  a  timing  uncertainty  °^T  =  15  ps,  corresponding  to 
150  chips,  due  to  position  uncertainty  (based  on  a  maximum 
position  uncertainty  of  +  5  n  mi) .  The  clock  update  procedure 
may  be  expected  to  take  a  maximum  of  about  20  sec  in  order  to 
collect  a  complete  frame  of  data  and  make  the  necessary 
computations.  The  time  to  search  the  t^me  uncertainty 

-60  /is  will  depend  on  the  receiver  in  question,  however  a 
maximum  search  time  of  about  43  seconds  would  be  required 
for  the  14  chip/sec  search  rate  quoted  for  the  results  of 
Figure  5 . 
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No  further  search  time  is  required  if  the  receiver  attempts 
to  acquire  the  additional  desired  signals  in  parallel. 

Thus  the  results  of  Figure  5  may  be  interpreted  as  follows. 
Hie  parallel  receiver  will  acquire  at  least  one  signal  in 
TacqC™*)  secon<*s  probability  .999.  Four  signals  may 
be  acquired  within  T  (max)  +  63  seconds  with  probability 
s  .9.  In  the  case  of  the  sequential  receiver,  the  first 
signal  is  acquired  within  TacCj(max)  with  probability  =  .9. 

A  maximum  of  about  twenty  seconds  is  then  required  to 
synchronize  the  clock,  and  a  maximum  of  43  seconds  is 
required  to  acquire  each  additional  signal  (with  probability 
^.9). 


The  results  of  Figure  6  may  be  interpreted  in  a  similar  manner 
except  that  about  100  seconds  is  required  to  acquire  additional 
signals  (due  to  the  6  chip /sec  search  rate) . 


In  the  interest  of  clarity,  we  will  for  the  moment  concern 
ourselves  only  with  the  parallel  receiver  which  acquires 
with  probability  ~.9  in  Tacq(max)  +  63  seconds  for  the  user 
with  velocity  uncertainty  <  5  MPH  and  in  Tac(j(max)  +  120 
seconds  for  the  user  with  velocity  uncertainty  160  MPH. 


5-6 


5.3 


Acquisition  Time  Conclusions 


Our  conclusions  may  be  summarized  as  follows: 

•  The  receiver  local  oscillator  is  the  limiting 
factor  in  performing  direct  acquisition. 

Direct  acquisition  appears  feasible  for  the  user 

with  velocity  uncertainty  <5  MPH  and  with  a  re- 

.0 

ceiver  local  with  overall  stability  <10  .  Direct 

acquisition  is  considerably  more  difficult  for 
the  user  with  velocity  uncertainty  as  great  as 
60  MPH.  In  this  case,  an  oscillator  with  stability 
on  the  order  of  10  ^  is  required  if  direct  acquisi¬ 
tion  is  to  be  practical  twenty-four  hours  after 
clock  synchronization. 

•  A  parallel  receiver  significantly  reduces  the 
time  required  to  acquire  four  signals. 

•  The  low  velocity  uncertainty  user  employing  a 
parallel  receiver  with  oscillator  stability 
<10"9  is  expected  to  acquire  all  four  signals 
within  about  300  seconds  with  probability  ^.9 
twenty-four  hours  after  clock  synchronization. 
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USER  OSCILLATORS 

The  characteristics  of  the  user  local  oscillator  are  an  important 
factor  in  the  ability  to  perform  direct  acquisition  The  quality 
of  the  user  oscillator  will  probably  vary  with  intended  applica¬ 
tion  since  better  performance  generally  implies  increased  size, 
weight  and  power  consumption.  For  the  purposes  of  this  investi¬ 
gation,  three  representative  oscillator  baselines  will  be  con¬ 
sidered: 

1.  A  low  power,  fast  warmup,  proportional  temperature 
control,  crystal  oscillator 

2.  A  high  quality,  double  proportional  control, 
crystal  oscillator 
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3.  An  atomic  standard  with  10  accuracy  (based 
on  the  User  Segment  direct  acquisition 
specification,  Paragraph  3.2.1.6  and  Table  II). 

A.l  Low  Power  Oscillator  (LPO)  Characteristics 

An  Army  sponsored  program  to  prove  feasibility  of  a  low  power, 
fast  warmup,  microcircuit  crystal  reference  oscillator  was 
conducted  by  Bendix  Corporation  and  is  reported  in  Reference  (3)  • 
Five  test  oscillators  were  built  under  this  program  with  varying 
degrees  of  success.  This  program  gives  an  indication  of  what 
one  may  expect  in  terms  of  oscillator  characteristics  for  a 
user  that  requires  direct  acquisition,  but  where  size,  weight, 
and  power  consumption  must  be  minimized.  The  program  goals 
and  best  results  achieved  are  summarized  in  Table  A-l.  Character¬ 
istics  of  a  commercially  available  oscillator  of  conventional 
design,  employing  a  proportional  control  oven,  are  also  included 
for  comparison.  A  baseline  oscillator  has  been  assumed  with 
characteristics  also  shown  in  Table  A-l.  The  baseline  in  most 

^H.H.  Greenhouse,  et  al,  "Fast  Warmup  Quartz  Reference  Oscillator, 
Final  Report,"  EGOM-0265-F,  Bendix  Communications  Div.,  June  1973. 
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instances  assumes  that  performance  will  be  at  least  as  good  as 
the  conventional  oscillator,  and  that  the  best  performance 
achieved  during  the  feasibility  program  will  be  repeated  on  a 
mass  production  basis  in  the  future. 

For  the  purposes  of  our  investigation,  let  us  assume  that  the 
user  oscillator  frequency  is  given  by 

fu  "  fo  +  *f  +ot 

where 

fQ  *  nominal  system  standard  frequency  (Hz) 

Af  °a  constant  user  frequency  bias  (Hz) 
a  *  frequency  aging  rate  (^fc)« 

The  frequency  bias,  Af,  will  be  considered  to  be  a  Gaussian 
random  variable  with  standard  deviation 

v  ■  £o  \R  *  ‘V2  +  ‘l2  ♦  *SIS2  ♦  *.2 

where 

■  temperature  stability 
m  voltage  stability 

Y  -  load  stability 

* -  short  term  stability 

STS 

■  setability 
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Substituting  the  parameters  assumed  for  the  baseline  oscillator: 


°Af  _ 

X" 


x  10'9)  +  (l0‘9)2  +(l0-9)2  +^2  x  10-10l 


=  6  x  10 
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^  -10  ^ 
+  10  y 

I 


Similarly,  a  will  be  considered  a  g.r.v.  with  1 6  value  equal  to 
the  aging  rate. 


The  baseline  low  power  oscillator  (LPO)  is  thus  characterised  by: 

\ 


=  6  x  10"9 


tt-  =  +  10*9/day. 

ZQ 


A. 2  Ultra-Stable  Crystal  Oscillator  (USO)  Performance 

Ttv~  temperature  stability  of  a  crystal  oscillator  can  be  improved 
by  approximately  an  order  of  magnitude,  at  the  cost  of  increased 
size  and  power  consumption,  if  a  double  proportional  controlled 
oven  is  employed.  The  characteristics  of  a  representative, 
commercially  available,  ultra-high  stability  oscillator  are 
shown  in  Table  A-2.  These  stability  characteristics  will  be 
assumed  as  a  baseline  for  a  sophisticated  user  who  does  not 
have  an  atomic  clock. 


A-4 


TABLE  A-2 


Vectron  Co-243 
Baseline  #2 

Size . 

4"  x  4“  x  9" 

Warmup  Power 

12w 

Operating  Power 

6W 

Short  Term  Stability 

<  10  ^/sec 

Frequency  Aging 

<  10_10/day 

Temperature  Stability 

+  10"10  0°C  to  50°C 

Voltage  Stability,  AV  =  5% 

+  2  x  10-11 

Load  Stability,  AR  -  10% 

+  10-11 

Frequency  Adjustment 

10-7 

Setability 

IQ’10 

Frequency  Recovery 

<  2  x  10'9 

Using  the  equations  of  Section  A.l,  the  baseline  ultra-stable 
crystal  oscillator  is  characterized  by 


x  10 


-10 


ta 


±  10'10/day 
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DIRECT  ACQUISITION  CODE  SEARCH  PERFORMANCE 

B.l  INTRODUCTION 

Direct  acquisition  of  the  P  component  of  the  Ll  signal  is 
being  considered  as  a  means  for  improving  the  ability  of  the 
GPS  user  to  operate  in  a  jamming  environment.  The  technique 
for  acquiring  the  P  code  is  to  step  a  correlator  in  1/2  chip 
increments  over  the  region  of  received  code  phase  uncertainty 
until  correlation  is  detected.  This  note  presents  analytical 
results  which  indicate  how  fast  the  code  can  be  searched  for 
a  given  set  of  performance  and  signal  parameters  for  the  non¬ 
coherent  correlator  of  Figure  B-l.  This  "in-lock"  detector, 
while  not  necessarily  optimum,  provides  good  performance 
consistent  with  a  minimum  of  complexity. 

The  IF  bandpass  filter  in  the  figure  is  centered  at  the 
nominal  frequency  of  the  input  signal,  and  has  a  bandwidth 
wide  enough  to  cover  the  frequency  uncertainty  of  the  input. 

The  output  of  the  IF  filter  is  fed  through  a  square- law 
device  whose  output  is  further  filtered  through  the  lowpass 
filter.  The  lowpass  filter  output  becomes  the  decision 
variable  based  on  which  the  code  acquisition  is  to  be  determined. 
That  is  "in- lock"  is  pronounced  when  the  decision  variable 
exceeds  a  certain  threshold  level \  otherwise  code  search  is 
continued. 


Figure  B-l  An  In-Lock  Detector 
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Referring  to  Figure  B-l,  let  the  input  signal  and  the  local 
oscillator  output  be 

r(t)  =  \/2P  s(t)  cos  (wQt  +  *>)  +  ^(t)  cos  (a>Qt  +  v>) 

-  r?2(t)  sin  (a>ot  +  *) 

q(t)  =  Vl  s(t  +  t)  cos  o>  t 

where 

s(t)  *  pseudonoise  (PN)  code  waveform  =  +  1 

s(t  +  €>  =  locally  regenerated  PN  code  waveform  with 
timing  error  of  M  <  a,  A  being  the  PN 
code  chip  time  interval, 

^(t),  2 (t)  *  two  independent  white  Gaussian  random 

processes  with  power  spectral  density 
of  NQ  watt /Hz. 

*  *  Uniformly  distributed  random  variable 
over  the  range  of  0  to  2n  radians. 

The  output  of  the  multiplier  is  given  by 

v(t)  *  VP  s(t)s(t  +«)  cos  (<tf^t  -  *) 

+  ±  s(t  +c)T?i(t)  cos  (cu.t  -*) 

y/2  1  1 

+  i  s(t  +€)T)  (t)  sin  (<*>.t  -*) 

V2  L 
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where  ok  =  -  a>0-  In  the  above  equation  we  have  ignored 

terms  involving  (*>  +  <o  )  because  they  fell  outside  of  the 

IF  filter  passband.  The  output  of  the  IF  filter  whose  band¬ 
width  is  very  much  smaller  than  the  PN  code  chip  rate,  is 

x(t)  =  Vp^L  “  cos  («^t  -*)  +  *7^)  cos  (^t  -*) 

+  sin  (wct  ) 

< 

when  gk  falls  inside  the  IF  filter  passband 

i  t 

n1<t)  cos  (Wjt  -¥>)  +  n2(t)  sin  (tOjt  -ifi) 
when  ok  falls  outside  the  IF  filter  passband 

f  f 

TK  (t)  and  V?(t)  are  independent  Gaussian  random  processes 

N 

with  power  spectral  density  of  watt/Hz. 

The  lowpass  component  of  the  square  law  device  output  is 

2  .  ,  p  /  |e|\2  lAt)  +  T)'2(t) 

?<■*>  ■  -2w|iowp.s8  ■  !  ■  V)  +  2—— 


+  /P  (l  -  ^)  Vt) 


when  falls  within  the  IF 
filter  passbaid 


|T,l^(t)  +  when  falls  out- 

1 - 2 - 

side  the  IF  filter  passband. 
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The  probability  density  of  y(t)  is  either  a  Rician  or  a 
Rayleigh,  depending  on  whether  the  signal  frequency  <o^  falls 
within  the  IF  filter  passband.  The  power  spectral  density  of 
y(t)  near  zero  frequency  is 


Sy(f)  = 


/p  N  Bt\  /  N  \  N 

'  (  j  +  -J-L  J  6  (f)+(p  +  Y-  Bj  J  when  <o±  falls 

within  the  IF  filter  passband 

,2 


(1) 


N  B 


N. 


when  oik  falls 


+\-f)  Bi 

outside  the  IF  filter  passband 


where  Bj  ■  IF  filter  bandwidth  in  Hz. 


The  decision  variable  z(t)  is  given  by 


(2) 


z(t)  -  y(t)*  h£(t) 

where  h^(t)  is  the  impulse  response  of  the  filter  and  stands 
for  convolution.  Ihe  probability  density  function  of  z  is  not 
readily  computed  because  y(t)  is  a  non-Gaussian  process. 

However,  if  the  bandwidth  of  the  lowpass  filter,  is 
much  smaller  than  that  of  the  IF  filter,  one  can  invoke  the 
central  limit  theorem  to  conclude  that  z(t)  can  be  approximated 
by  a  Gaussian  distribution.  That  is,  assuming  €  -  0 


pz(u) 


exp  - 


2°; 


2<u  "  mz) 
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where 


m  = 


P  .  "o"I 

2  2 

when  o), 

Vi 

when  co 

(3) 


°-Z2  = 


N 


N, 


p  +  "Y  Bi  )  -jr  when  co^  falls  within  the  IF  filter 


passband 


1  /  N  \  BjBjj  when  co^  falls  outside  the  IF  filter  passband 

1  \2/  1  (4) 


B  =  two-sided  noise  bandwidth  of  the  lowpass  filter 
N 


in  Hz 


(5) 


We  are  now  ready  to  compute  the  probabilities  of  false  alarm 
and  false  dismissal,  i.e.  probabilities  of  a  mistaken  "in- lock" 

when  in  fact  code  is  not  acquired  and  of  a  mistaken  "out-of-lock" 
when  code  is  acquired.  Letting  the  normalized  threshold  setting 
of  the  threshold  detector  be 


T  -  t 


/! 


where  t  being  actual  threshold  setting  in  volts,  then  it 
can  be  shown  that 


p .  *  probability  of  false  dismissal  -  Prob  (z  <  t  |  signal 


present) 


• 

r  i  "i 

1 

T  ' 

1  -  erf 

1  +  SNR  "  T 

,=i  erf  c 

1  +  — —  -  T 

1  SNR  1 

2 

~=  •  ^ 

7 

V«SNR/  SNR  BjJ 

L  V^  snr/snr  Bj  j 

(7) 
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P  =  probability  of  false  alarm  =  Prob(z>t  |  signal  not 
a 

present) 

=  j  erfc 

where 

SNR  *=  =  signal-to-noise  ratio  at:  the  IF  filter  output.  (9) 

woBI 

Pjj  *  (1  -  Pd)  and  Pa  have  been  computed  numerically  vs  SNR  for 
fixed  values  of  B^Bj  and  are  shown  in  Figures  B-2  through  B-5. 

In  the  analysis  presented  above  it  has  been  assumed  that  a  steady- 
state  condition  has  existed  in  the  system  prior  to  the  sampling 
of  the  decision  variable,  z(t),  and  a  decision  made.  That  is, 
the  dwell  time  TD  per  code  chip  is  at  least  several  times  longer 
than  the  time  constant,  Tc,  of  the  RC  lowpass  filter,  such  that 
the  signal  component  has  built  up  to  its  steady-state  value  at 
the  filter  output.  However,  for  the  problem  at  hand  we  are 
most  interested  in  minimizing  the  dwell  time  with  contraints 
on  the  values  of  Pa,  P^.B^Bj,  and  SNR.  Therefore,  we  have  to 
extend  our  analysis  to  include  the  case  where  the  system  is  in 
a  transient  state.  Again,  assuming  Bj^/Bj  «  1  such  that  the 
variable  z(t)  has  a  Gaussian  distribution,  the  mean  value  and 
variance  of  z  under  transient  conditions  can  be  shown  to  be 
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Figure 


Figure 
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N  B 

+  °~  when  falls  within  the 
IF  filter  passband 

when  o>^  falls  outside  the 
IF  filter  passband  (10) 


N  N 

+  -t^Bj  when  falls  within 

the  IF  filter  passband 

when  o>i  falls  outside  the 
IF  filter  passband  (11) 


The  respective  probabilities  of  false  dismissal  and  false  alarm 
are 


where  Tc  “  time  constant  of  the  RC  lowpass  filter  ■  •jgjj 
Figures  B-6  and  B-7  shows  (1  -  Pd)  vs  Td/Tc  with  SNR  and  Bj  as 
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Figure  B-7 


parameters.  In  both  figures,  Pa  is  assumed  to  be  0.005  and 

TD  =  0.1  and  0.05,  respectively.  It  is  seen  that  the  optimum 
value  of  Tp/T^  is  1.25  for  low  SNR  and  is  about  1.1  for 
values  of  SNR  such  that  PD  =  (1  -  Pd)  is  equal  to  or  greater 
than  0.9.  That  is,  for  maximum  probability  of  detection,  P^, 
the  dwell  time  should  be  1.1  times  the  time  constant  of  the 
lowpass  filter.  Figures  B-8  and  B-9  shows  P^  vs  (SNR)  for 
Bj  =  200  Hz  and  100  Hz,  respectively.  In  these  figures  Pa 
is  equal  to  .005  and  Tp/Tc  “  2  B^T^  is  set  at  1.125. 


PROBABILITY  OF  DETECTION 


CONCLUSION 


The  in-lock  detector  shown  in  Figure  1  has  been  analyzed.  It 
has  been  assumed  that  the  bandwidth  of  the  lowpass  filter  is 
much  narrower  thar  that  of  the  IF  filter  such  that  the  central 
limit  theorem  can  be  invoked  to  assert  that  the  decision  variable 
Z(t)  is  a  Gaussian  process.  The  probabilities  of  false  alarm 
and  false  dismissal  (Pa  and  Pd)  have  been  derived  for  both 
transient  and  steady-state  conditions  of  the  system.  Numerical 
results  of  (1  -  Pd)  vs  several  different  system  parameters 
have  been  obtained.  It  has  been  shown  that  the  optimum  value 
of  T  /T  (dwell  time/time  constant  of  the  lowpass  filter) 
for  maximum  value  of  (1  -  Pd)  is  approximately  1.1  for  the 
range  of  (SNR)  being  considered. 
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GPS  SPECIAL  STUDIES  AND  ENGINEERING  PROGRAM 


TASK  III  -  MANPACK  STUDY 


1 .0  INTRODUCTION 

This  report  presents  the  results  of  a  study  concerning  the 
proposed  Army  manpack  GPS  user  to  be  operational  in  the  early 
1980's.  The  objectives  of  the  study  were  to: 

1)  Provide  the  Army  an  insight  to  the 
compatibility  of  the  user  established 
performance  characteristics  as  set 
forth  in  the  Army  Mission  Requirements 
Document,  YEN-74-98,  with  the  imple¬ 
mentation  of  the  GPS  manpack. 

2)  provide  tradeoff  information  to 
assess  the  impact  of  any  proposed 
relaxation  of  goals 

3)  Provide  a  recommended  approach  to  the 
manpack  effort. 

The  report  contains  three  major  sections.  Section  2  deals  with 
the  manpack  overall  equipment  configuration  and  performance 
recommendations.  Section  3  discusses  microprocessor  hardware 
and  Section  4  deals  with  computer  functions  and  software  recom¬ 
mendations  . 
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1.1 


Study  Conclusion  Summary 


The  study  conclusions  may  be  summarized  as  follows: 

1)  The  manpack  goals  (Table  1-1)  and 
scenario  do  not  lead  to  any  major 
receiver  hardware  simplification  when 
compared  to  other  users  (equipment 
interfaces  are  simplified  due  to  the 
absence  of  external  sensors) . 

2)  A  manpack  set  which  satisfies  the  mission 
requirements  of  a  majority  of  users 

as  stated  in  the  nArmy  Mission  Require¬ 
ments  Document”  can  probably  be  achieved 
within  the  unit  size,  weight,  and  cost 
goals  (Table  1-1)  given  extensive  LSI 
development . 


Table  1-1  Manpack  Development  Goals  Summary 


Horizontal  Accuracy  (CEP) 
Vertical  Accuracy  (CEP) 
User  Motion 
Interference  Rejection 
Mission  Duration 
Weight 
Volume 
Unit  Cost 


1.5  to  50  m 
1.1  to  50  m 
0  to  50  MPH 
High  AJ  Desirable 
2  to  48  Hours 
8  to  12  lbs. 

500  to  1000  cu.  m. 
$15,000 
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3)  A  "full  performance”  dual  manpack/ 
vehicular  unit  will  probably  not  meet 
the  size  weight,  and  cost  goals  set 
forth  in  Table  1-1,  but  appears  feasible 
if  some  relaxation  of  goals  is  acceptable. 

4)  A  low-cost,  low-power  computer  capable  of 
performing  the  required  navigation  computations 
for  the  low  performance  user  is  feasible  using 
available  microprocessors;  however  the  high 
performance  user  will  require  either  higher 
speed,  hence  higher  power,  bi -polar 
technology  or  augmentation  of  existing 
microprocessors  with  external  floating  point 
hardware.  The  high  accuracy  user  can  be 
accomnodated  at  the  expense  of  computation 
speed  through  the  user  of  multiple-precision 
arithmetic . 

5)  Because  of  the  shorter  word  length  associated 
with  the  microprocessor  the  navigation  filter 
must  be  implemented  as  a  "square -root”  filter 
to  improve  accuracy  and  guarantee  stability. 


S t udy  Recommendation  Summary 


The  study  recommendations  may  be  summarized  as  follows: 


1)  It  is  recommended  that  the  Army 

consider  developing  separate  low-cost 
manpack  (LCM)  and  high-performance 
manpack/vehicle  (HPM)  user  equip¬ 
ments.  The  anticipated  equipment 
characteristics  are  summarized  in 
Table  1-2. 
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2)  "All-digital"  receiver  implementation 
should  be  pursued  with  a  maximum  of 
operations  performed  by  a  micro -computer . 

3)  LSI  development  of  common  functional 
elements  should  begin  as  quickly  as 
possible. 

4)  Thick  film  development  of  analog 
circuitry  should  be  investigated. 

5)  Additional  investigation  into  hardware/ 
software  trade-offs  in  receiver  control, 

code  generation,  data  demodulation  and  floating 
point  computation  to  determine  the  optimum 
division  of  functions  between  the  processor 
and  dedicated  hardware. 

6)  Consider  further  the  feasibility  of  "add-on” 
hardware  floating-point  module  to  increase 
computation  speed  of  the  microprocessor  for 
use  in  high-performance  applications. 

7)  Recommendations  -  The  navigation  filter  for 
the  manpack  user  should  be  a  low-order  Kalman 
filter  implemented  in  square-root  form  (see 
Section  4.3.1  and  Appendix  C) . 


Table  1-2 


MANPACK  CHARACTERISTICS  SUMMARY 


High  Performance 
Manpack/Vehicle  Set 
Characteristics  Summary 

Low-Cost 

Manpack  Set 

Characteristics  Summary 

•  Dual  raanpack/vehicle  use. 

•  Manpack  use  only  (0-5  mph). 

•  10  m  (CEP)  position 

•  20  m  horizontal  position 

accuracy. 

accuracy  (CEP). 

•  Maintain  accuracy  in  the 

•  10  m  vertical  position 

presence  of  40  dB  of 

accuracy  in  "high  vertical 

jamming. 

resolution*'  mode. 

#  TTFF  on  the  order  of  1 

•  Maintain  accuracy  in  the 

minute. 

presence  of  30  dB  of  jamming. 

•  12  hour  continuous  operation 

•  TTFF  on  the  order  of  4 

on  internal  batteries. 

minutes. 

•  23  lbs.  total  weight. 

•  12  hours  of  continuous 

operation. 

•  12  lbs  total  weight. 

Anticipated  Cost 

Anticipated  Cost 

$25,000  to  $35,000 

$10,000  to  $15,000 

(1000  production  units) 

(1000  Production  units) 
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SECTION  2 


2.0  MANPACK  CONFIGURATION 


2.1  Manpack  Development  Goals 

It  is  the  Army's  objective  to  develop  a  small,  lightweight 
user  equipment  for  manpack  as  well  as  wheel  and  track  vehicle 
use.  The  manpack  user  equipment  (MUE)  is  intended  to  satisfy 
needs  for  a  wide  variety  of  missions  (see  Ref.  1  )  and  would 
preferably  require  that  only  one  equipment  be  developed  in 
order  to  achieve  lowest  acquisition  and  maintenance  costs. 
Present  MUE  goals  are: 

Weight:  8-12  pounds 
Volume:  500  -  1000  cubic  inches 
Mission  Duration:  48  hours 
Unit  Cost:  $15,000 


Required  performance  is  mission-dependent  and  varies  considerably. 
The  range  of  performance  is  summarized  below. 


Horizontal  Accuracy  (CEP) 
Vertical  Accuracy  (CEP) 
User  Motion 

Interference  Rejection 
Mission  Duration 


1.5  to  50  m 
1.1  to  50  m 
0  to  50  mph 
High  AJ  desirable 
2  to  48  hours 


An  example  manpack  set  is  shown  in  Figure  1. 

It  is  obviously  difficult  to  accurately  predict  the  ability 
to  meet  these  goals  in  five  years  time  in  the  rapidly  changing 
technology  and  economic  situations  of  today.  Unit  end  cost 
is  especially  difficult  to  predict  since  a  low  per  unit  cost 
may  be  achieved  only  through  high  development  costs.  For 
example,  the  size,  weight,  and  cost  reductions  exemplified. 

Ref  (1)  MArmy  Mission  Requirements  Document ,M  YEN-74-98. 
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by  today's  microprocessors  were  realized  through  extensive 
development  programs  justified  by  the  vast  potential  market 
and  the  "universal"  character  of  the  computer. 

In  view  of  the  uncertainties  involved,  we  have  tried  to  be 
relatively  conservative  in  our  estimates  of  what  can  be 
accomplished.  Further,  the  recommendations  are  intended  to 
provide  direction  for  developments  which  will  lead  to  the 
most  satisfactory  solutions  regardless  of  absolutes. 

2.2  Recommended  User  Classes 

Present  circumstances  do  not  justify  the  assumption  that  one 
manpack /vehicle  UE  which  satisfies  all  the  mLssion  requirements 
can  be  developed  by  1980  within  the  size,  weight  and  cost 
goals  stated  above  (for  a  total  production  of ~5000  units). 

Of  these,  the  cost  goal  would  seem  to  be  the  most  difficult 
to  achieve  for  a  unit  with  the  required  performance  and  close 
to  the  desired  size  and  weight. 

The  ultimate  manpack  goals,  as  presently  stated,  are  such 
that  required  complexity  appears  to  be  approximately  the  same 
as  the  set  "X"  or  as  a  minimum,  the  set  "Y"  presently  under 
development.  The  low  user  dynamics  (0  -  50  mph)  do  not  seem 
to  allow  any  major  simplification  other  than  that  external 
sensors  are  not  necessary.  Thus  the  manpack  receiver  requires 
a  microminiaturized,  low  power  version  of  a  complex  receiver. 
Such  a  receiver  will  probably  not  come  at  low  cost,  even  in 
relatively  large  quantity.  We  would  expect  the  price  to  be 
somewhere  between  $25  and  $35  thousand  dollars  per  unit  in 
quantities  of  1000  units. 

The  end  price  depends,  of  course,  on  the  number  of  units  pur¬ 
chased  from  a  single  supplier  at  one  time  as  well  as  total 
expected  production.  A  buy  of  100  units,  which  is  probably  a 
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more  realistic  quantity  in  the  early  stages  of  system  opera¬ 
tion,  would  result  in  considerably  higher  cost,  probably  in 
the  range  of  $35K  to  $50K  (excluding  any  development  costs). 

Many  of  the  missions  for  which  the  manpack  is  intended  have 
requirements  which  could  be  met  with  a  UE  of  significantly 
reduced  complexity  and  cost.  In  particular,  it  appears  that 
Ll  C -signal  only  operation  using  a  single  channel,  slow 
sequencing  receiver  could  provide  position  accuracy  on  the 
order  of  20  meters,  for  a  user  with  a  0  -  5  mph  velocity, 
under  most  conditions.  The  estimated  price  of  such  a  receiver 
is  on  the  order  of  $10  to  $15  thousand  dollars  in  quantities 
of  1000  ($15  to  $20  thousand  dollar  ,  in  lots  of  100  excluding 
nonrecurring  costs).  The  major  disadvantage  is  the  reduction 
in  AJ  performance.  Thus  it  is  recommended  that  the  Army  con¬ 
sider  development  of  two  manpack/vehicular  UE  as  follows: 

•  A  full  capability  equipment  for  dual 
manpack/vehicular  use 

•  A  Ll  C-signal  only  receiver  for 
manpack  use 

A  more  complete  discussion  of  the  development  recommendations 
and  problems  is  given  below. 

2.3  High  Performance  Manpack 

The  High  Performance  Manpack  (HPM)  is  intended  to  meet  stringent 
mission  requirements  with  design  goals  and  features: 

•  Dual  manpack /vehicle  use 

•  10  m  position  accuracy  (CEP) 

•  Maintain  accuracy  in  the  presence  of 
40  dB  of  jamming 

•  TTFF  on  the  order  of  1  minute 

•  12  hr  continuous  operation  on  internal  batteries 

•  23  lbs  total  weight 
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The  above  goals  imply 


•  L1/L2  operation 

•  P-signal  operation 

The  complexity  of  such  an  equipment  is  assumed  to  be  approxi¬ 
mately  that  of  the  "X-Set"  equipment  presently  under  develop¬ 
ment.  The  somewhat  reduced  dynamic  environment  does  not  seem 
to  allow  any  major  simplification  other  than  that  external1 
sensors  are  probably  not  required. 

The  receiver  configuration  would  ideally  be  a  fully  continuous 
4-channel  receiver,  however  this  is  not  a  likely  choice  unless 
new  developments  in  the  code  tracking  loop  circuitry  bring 
about  a  reduction  in  size  and  cost.  Another  potential  con¬ 
figuration  provides  continuous  carrier  tracking  but  employs  a 
single  time -sequenced  code  loop.  This  is  the  technique  being 
implemented  for  the  present  MX”  set.  The  third  choice  of  a 
single  code/carrier  channel  time-sequenced  over  the  various 
signals  has  considerably  reduced  AJ  performance  and  is  not 
recommended  for  the  HPM, 


2.4  Low-Cost  Manpack 

The  recommended  low-cost  manpack  (LCM)  is  basically  a  C-signal 
only,  slow  sequencing  receiver,  designed  to  provide  moderate 
performance,  adequate  for  many  missions,  with  a  minimumt.of: 
equipment  size  and  costv 

The  use  of  the  C-^ignal  simplifies  both  circuitry  and  opera¬ 
tional  difficulties.  This  is  gained  at  the  expense  of  accuracy 
and  AJ  performance.  The  ranging  accuracy  is  not  greatly 
reduced  from  the  ? -signal  for  a  low  dynamic  user.  Range 
accuracies  of  10  feet  should  be  easily  achieved  for  the  man- 
pack  user.  Ionospheric  errors  are  a  major  source  of  errors 
if  the  C-signal  is  transmitted  on  Ll  only.  Simple  modeling 
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can  limit  the  error  to  about  15  feet  under  most  circumstances. 

It  is  possible  that  this  error  can  be  further  reduced  by 
inputting  predicted  delays  for  missions  of  short  duration, 
more  complicated  models,  or  meteorlogical  data  received  through 
an  auxiliary  link.  Figure  2  shows  a  typical  situation. 

The  Army  missions  requirements  (Ref.  1  )  indicate  that  the 
LCM  would  be  especially  useful  if  a  vertical  position  accuracy 
on  the  order  of  10  m  can  be  achieved.  Thus  a  "high  vertical 
resolution"  mode  is  suggested  in  which  satellite  selection 
and  position  computation  are  optimized  in  the  vertical  direction. 

The  LPM  is  intended  to  have  the  design  goals 

•  20  m  horizontal  position  accuracy  (CEP) 

•  10  m  vertical  position  accuracy  in 
"high  vertical  resolution"  mode 

•  Maintain  accuracy  in  the  presence 
of  30  dB  of  jamming 

•  TTFF  on  the  order  of  4  minutes 

•  Manpack  use  only  (0-5  mph) 

•  12  hrs  of  continuous  operation 

•  12  lbs  total  weight 

It  is  anticipated  that  the  LCM  would  be  implemented  as  a  slow 
sequencing  receiver  (similar  to  the  "Z"  set  presently  under 
development).  Some  of  the  advantages  and  disadvantages  of 
the  LCM  are  listed  in  Table  2-1.  A  skeleton  specification 
for  the  LCM  is  included  as  Appendix  A. 
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Figure 


Auxiliary  Data  Link  to  Reduce 
LCM  Ionospheric  Errors 


Table  2-1  LCM  Features 


FEATURES 

•  Single  frequency  operation 

•  No  P-Code  generators  or  acquisition 
circuitry  required 

•  No  secure  code  key  or  "HOW"  required 

•  Rapid  reacquisition 

•  Single  channel  code  and  carrier 
tracking  circuitry 

•  Fully  automatic  and  "autonon.ous"  operation 
DISADVANTAGES 

•  Sequential  operations  increase  TTFF 

•  Interference  resistance  reduced  by  about  10  dB 

•  Resulting  accuracy  of  20  to  30  meters 
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Comparison  of  ar>d  LCM 


The  recommendation  to  develop  two  tnanpack  units  is  predicated 
on  the  assumptions  that 


A  high  performance  receiver  will  cost 
significantly  more  than  $15K  even  in 
quantities  of  1000  units 

The  need  to  accomplish  missions  requiring 
high  performance  is  such  that  the  cost  of 
the  equipment  is  justified 

A  much  simpler  and  lower  cost  unit  can  be 
developed  which  will  meet  the  needs  of  a 
1  arop  number  of  users 


The  LCM  savings  are  gained  through  reduced  accuracy,  near 
stationary  user  (-unpack  only),  less  »J  protection,  and 
longer  TTFF  as  discussed  in  the  preceding  sections. 

The  possible  simplifications  may  be  summarized  as  follows.  The 
reduced  accuracy  and  AJ  performance  allows 


ingle  frequency  operation  (increased 
onospheric  error) 

-signal  tracking  only  (lower  ranging 
°  -,,,1  at  orotection) 


n/UirOPU 


simpler  software  and  CPU  (lower 
computational  accuracy) 


The  near  stationary  user  (  5  mph)  and  longer  TTFF  leads  to  the 
possibility  of 


nele  channel  slow  sequencing  operation 
JL.,#-  comolex  control  hardware 


and  software 
no  external  sensors 


soma  Insight  Into  the  computation  complexity  of  the  too  agulp- 
I  Vf  ghihdl  'tom  Figure  3  .  The  overall  block  diagram 

represents  the  HPM  while  the  shade  element.  Indie.te  those 

required  for  the  L£M. 
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FIGURE  3 

CONFIGURATION  FOR  COMPARING  HPM  AND  LCM  COMPLEXITY  DEMOD 
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Power  Dissipation 


Power  dissipation  is  one  of  the  most  critical  design  areas, 
since  the  battery  pack  is  a  major  contributor  to  the  UE 
size  and  weight.  Further,  heat  dissipation  in  the  small 
volume  of  the  manpack  can  have  a  critical  effect  on  component 
temperature.  For  example,  the  "Y"  set  presently  under  develop¬ 
ment  is  estimated  to  require  543  watts  of  power  including  45 
watts  for  six  cooling  fans . 

Lithium  batteries  presently  under  development  represent  a 
major  breakthrough  in  battery  technology,  providing  about 
an  order  of  magnitude  improvement  in  energy  storage  per 
pound  of  battery  weight  over  other  rechargeable  batteries 
(see  Appendix  B).  Assuming  a  battery  capacity  of  150  w-hr/lb 
for  estimation  purposes,  a  pound  of  battery  is  required  for 
every  12.5  watts  of  equipment  power  for  a  12  hour  mission. 
Three  pounds  of  batteries  could  thus  provide  37.5  watts  of 
continuous  power.  Low  power  logic  technology  is  rapidly 
improving  which  will  help  in  meeting  this  design  requirement. 

A  major  reduction  in  power  consumption  is  achievable  through 
judicious  use  of  power  switching  of  components  that  do  not 
require  continuous  operation.  There  is,  at  present,  about 
a  4  to  1  difference  in  power  consumption  for  digital  logic 
with  maximum  operating  speeds  of  around  10  Mbps  and  1  Mbps. 
Thus  there  is  apparent  advantage  to  C-signal  only  operation 
in  addition  to  reduced  complexity. 

The  HPM  is  expected  to  require  about  45W  of  battery  power. 

Note  that  this  is  approximately  an  order  of  magnitude  less 
than  the  present  Y-set  development  model.  This  reduction  in 
power  is  hoped  to  be  gained  through  careful  attention  to 
circuit  configuration  for  minimum  power  dissipation,  use  of 
low  power  logic  and  components,  and  power  switching. 

The  LCM  is  expected  to  require  about  25W  of  battery  power. 

The  reduction  over  the  HPM  is  due  to  less  complicated 
circuitry  as  well  as  lower  clock  rates. 
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Microminiaturization 


The  size  goal  for  a  manpack  receiver  is  to  achieve  a  volume 
of  less  than  1000  cu.  in.  including  the  battery  pack.  This 
corresponds  to  something  greater  than  a  5  to  1  reduction  in 
volume  of  the  Y  set  presently  under  development.  A  reduction 
of  this  magnitude  does  not  seem  entirely  out  of  the  question 
using  existing  techniques,  as  may  be  inferred  from  the  reduction 
in  computer  size  that  is  realized  with  the  microprocessors 
now  available.  It  must  be  remembered  however  that  the  extensive 
LSI  developments  required  were  undertaken  because  of  the 
'’universal"  nature,  and  therefore  large  potential  market,  of 
the  device.  The  state-of-the-art  and  availability  of  micro¬ 
processors  is  such  that  they  should  be  employed  to  the  maxi¬ 
mum  extent  possible  in  place  of  special  purpose  LSI  components. 
Where  this  is  not  possible,  due  to  processing  speed,  it  would 
seem  most  cost-effective  to  develop  LSI  devices  to  be  made 
available  to  all  UE  manufacturers,  e.g.  single  LSI  code 
generator  chip  would  be  developed  as  opposed  to  each  manu¬ 
facturer  attempting  to  develop  his  own.  In  order  to  take 
advantage  of  the  rapid  progress  in  digital  circuit  technology, 
a  maximum  of  the  receiver  functions  should  be  implemented 
digitally. 

The  analog  circuitry  is  probably  best  handled  with  thick 
film  techniques,  with  the  possible  exception  of  the  L-band 
circuitry,  where  thin  film  may  be  more  advantageous.  While 
it  appears  that  thick  film  implementation  of  the  IF  and  low 
frequency  analog  circuitry  is  practical,  experimental 
development  work  is  required  to  substantiate  this  conclusion. 
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Weight 


Equipment  weight  depends  largely  on  the  extent  to  which  LSI 
and  microminiaturization  analog  circuitry  is  employed,  as 
well  as  the  module  packaging  and  functional  integration 
employed,  i.e.  if  several  functions  are  integrated  in  one  chip, 
and  the  number  of  separate  module  covers  and  interconnectcrs 
are  minimized,  the  weight  will  be  reduced  over  a  receiver 
employing  a  lower  level  of  chip  and  module  integration.  Power 
consumption  also  affects  weight  since  the  battery  weight  is 
considered  as  part  of  the  equipment  weight. 

The  HPM  is  expected  to  weigh  about  23  lbs.  including  4  lbs. 
of  batteries,  while  a  total  weight  of  12  lbs.  seems  reasonable 
for  the  LCM  (including  2  lbs.  of  batteries). 


2.9  Receiver  Hardware 

The  manpack  effort  is  primarily  one  of  microminiaturization. 
In  order  to  be  successful,  the  Army  program  should: 

•  Exploit  microprocessors  to  the 
fullest 

•  Start  LSI  development  of  certain 
common  functional  blocks 

•  Start  thick  film  development  of 
analog  circuitry. 

The  basic  block  diagram  for  the  manpack  user  equipments  (UE) 
is  shown  in  Figure  4.  A  more  detailed  block  diagram  of 
a  single  correlator,  carrier  tracking  loop,  noncoherent  code 
tracking  loop,  and  data  demodulator  is  shown  in  Figure  5. 
There  are,  of  course,  many  versions  of  this  circuitry 
and  this  particular  one  is  used  only  as  an  example. 
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Simplified  Manpack  Set  Block  Diagram 
Figure  4 
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This  circuitry  and  the  associated  control  logic  represents 
most  of  the  complex  circuitry  in  the  receiver  and  is  a  candi¬ 
date  for  nearly  all-digital  implementation.  A  block  diagram 
of  the  digital  version  is  shown  in  Figure  6.  The  dotted 
lines  represent  functions  that  can  probably  be  performed  by 
a  microprocessor  while  the  crosshatched  blocks  represent 
digital  LSI. 

It  is  desireable  to  perform  as  many  operations  as  possible 
with  microprocessors,  due  to  the  availability  and  cost.  Some 
functions  must  be  done  in  LSI  firmware,  however,  due  to  pro¬ 
cessing  speed  limitations. 

The  sampling  rate  and  the  number  of  bits  used  in  A/D  conversion 
must  be  chosen  such  that  the  digital  processing  introduces 
negligible  degradation.  These  parameters  require  careful 
selection  in  the  design  process.  The  data  detection  filter 
must  process  samples  at  about  10  times  the  data  rate  (i.e.  500 
samples/sec)  in  order  to  introduce  less  than  0.5  dB  of 
degradation  (Ref.  2  ).  Experience  has  shown  that  4  bits  of 

quantization  (16  levels)  is  usually  adequate  to  ensure  negli¬ 
gible  degradation  in  BER  due  to  quantization  noise.  Four-bit 
(16  level)  quantization  is  adequate.  The  carrier  tracking  loop 
filter  can  operate  at  a  somewhat  lower  rate,  on  the  order  of 
100  to  200  samples/second  for  loop  bandwidth  from  20  to  40  Hz. 
Three  bit  quantization  is  recommended  although  fewer  bits 
could  be  used.  The  code  tracking  loop  need  only  operate  at 
around  20  to  50  samples/sec  due  to  the  very  low  bandwidth 
(<10  Hz)  involved.  All  of  the  above  rates  are  within  micro¬ 
processor  capabilities.  Some  of  the  other  functions,  such  as 
code  generation  and  the  digital  number -controlled  oscillator 
(NCO)  involve  high  clock  rates  and  probably  will  require 
dedicated  LSI  hardware. 

Ref  (2)  F.D.  Natali,  "All-Digital  Coherent  Demodulator  Techniques 
ITC  Proc.,  Los  Angeles,  California  1972. 
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Digital  Implementation  of  Carrier 
and  Code  Tracking  Loops 


due  to  the  very  low  bandwidth  (<10  Hz)  involved.  All  of  the 
above  rates  are  within  microprocessor  capabilities.  Some 
of  the  other  functions,  such  as  code  generation  and  the  digital 
number-controlled  oscillator  (NCO)  involve  high  clock  rates 
and  probably  will  require  dedicated  LSI  hardware. 

The  correlator,  as  shown,  is  implemented  in  an  analog 
fashion.  This  approach  is  recommended  over  a  digital  imple¬ 
mentation  with  present  technology  since  digital  correlation 
involves  high  clock  rates  as  well  as  considerable  complexity, 
and  there  does  not  appear  to  be  any  particular  advantage  to 
be  gained  since  most  of  the  analog  circuitry  is  passive. 

The  HPM  would  ideally  incorporate  four  signal  tracking  and 
data  demodulation  channels.  Microprocessor  architecture  is 
such  that  similar  functions  can  be  performed  for  all  four 
channels  by  a  single  processor  if  adequate  speed  is  available. 
As  noted  above,  a  single  channel  requires  relatively  simple 
processing  of  about  700  samples /second  and  four  channels 
would  thus  require  processing  about  2800  samples /second.  The 
sampling  requirements  are  summarized  in  Table  2-2. 

The  functions  denoted  as  LSI  hardware  require  higher  speeds 
and  are  not  expected  to  be  accomodated  by  computer-type 
architecture.  The  question  then  arises  as  to  whether  size 
or  cost  savings  can  be  achieved  by  time-multiplexing  certain 
high  speed  LSI  elements.  The  code  generators  do  not  seem  to 
be  significantly  simplified  by  time-sharing  common  circuit 
elements  and  independent  implementation  seems  the  best  choice. 
The  NCOt  on  the  other  hand,  may  benefit  from  time-sharing 
of  certain  elements,  depending  on  the  technique  employed.  A 
possible  block  diagram  of  a  four-channel  receiver  employing 
time  multiplexed  hardware  issshown  in  Figure  7. 
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Function 

Bandwidth 

Recommended 

Sample  Rate 

Quanti¬ 

zation 

Data  Detection 

Rd  =  50  Bps 

^500  samples /sec 

4  Bits 

Carrier  Tracking 
Loop 

Bl  <  40  Hz 

<200  samples/sec 

3  Bits 

Code  Tracking 
Loop 

Bl  <10  Hz 

<50  samples/sec 

3  Bits 
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3.0  MICROPROCESSORS  AND  MICROCOMPUTERS 

3.1  Introduction 

The  size,  power  and  weight  requirements  for  the  manpack  unit 
influence  to  a  great  extent  the  choice  of  computer.  Because 
of  these  factors,  use  of  existing  aerospace  computers  in 
manpack  operations  would  greatly  compromise  the  design. 
Fortunately  certain  of  the  microprocessors  presently  being 
developed  offer  smaller  size,  reduced  cost  and  much  lower 
power  consumption. 

In  the  present  baseline  design  the  computer  is  responsible 
for  both  calculation  and  receiver  control,  although  separate 
circuits  for  these  functions  are  being  considered.  At  the 
present  time  a  16  bit  binary  machine  with  floating  point 
capability  appears  to  provide  acceptable  accuracy  (6  to  7 
digits)  as  long  as  procedures  which  are  tolerant  of  round-off 
errors  are  employed,  particularly  in  the  Kalman  filter  routine. 

For  higher  precision  one  alternative  to  going  to  24  or  32 
bit  machines  or  to  double  precision  software  on  a  16  bit 
machine  is  to  utilize  BCD  arithmetic,  which  allows  essen¬ 
tially  unlimited  computational  accuracy  given  enough  time. 
Trading  increased  computation  time  for  reduced  hardware 
requirements  appears  to  be  one  area  that  can  yield  consider¬ 
able  benefits  for  the  manpack  design. 

There  are  a  number  of  microprocessors  either  available  or 
soon  to  be  available  which  would  have  generally  suitable 
characteristics,  including  the  Teledyne  TDY52A  and  TDY52B, 
the  Intel  8080,  the  Intersil  6100  (PDP-8  code  and  CMOS  con¬ 
struction),  the  RCA  COSMAC  (also  CMOS)  and  the  National 
Semiconductor  IMP-16.  The  major  difficulties  with  using 
these  microprocessors  is  the  lack  of  higher  level  languages 
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(except  for  the  Teledyne,  Intel  and  Intersil  machines)  and 
their  special  interfacing  requirements.  One  feature  that 
must  be  carefully  examined  is  the  ability  to  mix  RAM  and 
ROM  memory  and  the  capability  of  addressing  sufficient 
memory  (at  least  16K)  to  hold  all  of  the  required  programs. 
Another  feature  that  is  particularly  important  for  manpack 
use  is  external  interrupt  capability  to  allow  the  receiver 
to  signal  the  computer  that  receiver  functions  are  complete. 
Other  characteristics  that  must  be  carefully  examined  are 
the  size  and  power  consumption  of  the  total  package  with  all 
required  memory  and  interface  chips,  along  with  the  other 
considerations  of  speed,  second  sourcing  and  price. 

One  of  the  interesting  developments  that  may  be  of  great 
significance  in  the  early  1980  time-frame  is  the  recent 
introduction  of  computer  "slices."  These  chips  are  usually 
either  2  or  4  bits  of  a  complete  CPU  which  may  be  stacked 
together  to  form  a  computer  with  any  desired  word  length, 
freeing  the  designer  to  select  the  particular  word  length 
which  suits  the  problems  and  reducing  the  complexity  of 
software  required  to  maintain  suitable  accuracy  with  short 
word  length  computers. 

3.2  Review  of  Present  Microprocessor  Technology 

At  the  present  time  the  number  of  microprocessor  types  is 
rapidly  increasing  with  at  least  26  different  models  either 
presently  available  or  available  in  the  near  future.  Intel 
Corporation  has  five  models  -  two  of  them  are  4-bit  micro¬ 
processors  to  8-bit  microprocessors  and  a  bi-polar  slice 
kit  which  can  be  used  to  make  computers  up  to  32  bits  wide. 
National  Semiconductor  has  five  models  -  three  of  them  in 
their  IMP  line,  one  called  GPC/P,  and  a  recently  introduced 
single  chip  16-bit  model  called  PACE.  Rockwell  has  both  a 
4-bit  and  an  8-bit  chip  set.  Both  RCA  and  Intersil  have 
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complementary  MDS  (CMOS)  chips  which  are  attractive  because 
of  the  very  complete,  interfacing  chips  that  are  available 
for  use  with  the  basic  CPU. 

In  addition  to  these  manufacturers,  Fairchild,  Signetics, 

Electronic  Arrays,  Micro  Systems  International  and  Toshiba 

also  produce  an  n-channel  device,  while  Mostek,  Burroughs, 

and  Fairchild  offer  p-channel  systems.  Monolithic  Memories, 

Raytheon  and  Transitron  build  4-bit  microprocessor  slices 

using  bi-polar  technology.  Texas  Instruments  is  planning 

to  introduce  a  4-bit  slice  using  current -inject ion  logic 
o 

(I  L)  which  has  the  unique  property  that  power  consumption 
can  be  traded  for  increased  speed  simply  by  increasing  the 
supply  current.  Actron  Industries  is  presently  producing 
a  control  unit  and  an  arithmetic  unit  as  costom  chips  for 
a  16-bit  machine  which  can  be  expanded  in  eight  bit  slices 
up  to  80  bits  for  use  as  a  standardized  avionics  module  for 
the  Air  Force.  Honeywell  is  developing  a  complete  micro¬ 
computer,  the  MOD/LSI-2 ,  for  use  in  small  missiles,  which 
may  be  suitable  for  manpack  use  and  because  of  its  large 
volume  use  in  missiles  may  be  reasonably  priced  (however, 
very  little  data  on  this  device  is  available). 

3.3  System  Considerations  in  Microprocessor  Selection 

At  the  present  time,  most  of  the  single  chip  microprocessors 
are  considerably  slower  and  less  flexible  than  existing 
minicomputers.  Comparison  of  the  execution  time  for  several 
microcomputers  and  minicomputers  show  that  the  execution 
time  can  be  between  about  2  and  4  times  slower  in  the  case 
of  the  microcomputer.  The  primary  reason  for  this  is  that 
the  microprocessors  typically  are  based  on  metal  oxide 
semi-conductor  technology  while  the  minicomputers  are  built 
using  bi-polar  devices.  Another  important  reason  for  the 
difference  in  speed  is  that  most  of  the  microprocessors  use 
a  smaller  number  of  bits  -  typically  either  4  or  8  bits. 

One  way  of  comparing  the  speed  of  microprocessors  and  mini- 
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computers  is  to  measure  the  execution  time  for  typical 
functions  which  would  be  performed  in  the  system.  Compari¬ 
son  of  execution  speeds  of  an  8  bit  Intel  8080  microprocessor 
which  is  representative  of  types  presently  available  on 
the  market  and  an  Inter-Data  50  minicomputer  which  is  a 
typical  16-bit  bi-polar  minicomputer,  shows  that  for  a 
typical  error  check  operation,  such  as  for  a  cyclic  redun¬ 
dancy  check  code,  the  Intel  8080  takes  somewhere  between 
250  and  320  microseconds  to  execute  whereas  the  Inter-Data 
50  will  execute  that  same  set  of  instructions  in  12  to  15 
microseconds.  The  reason  for  this  discrepancy  in  execution 
times  can  be  traced  to  the  relative  slowness  in  executing 
individual  instructions  of  the  microprocessor.  Again,  com¬ 
paring  the  Intel  8080  and  the  Inter-Data  50,  subtraction  or 
addition  takes  about  2  microseconds  for  the  8080  and  only  1 
microsecond  for  the  Inter-Data,  register  loads  take  about 
3.5  and  1  microsecond  respectively,  and  memory  reference 
instructions  require  about  3.5  microseconds  and  3.25  micro¬ 
seconds  respectively.  For  these  relatively  simple  instruc¬ 
tions  the  microprocessor  and  the  minicomputer  compare  quite 
favorably  with  only  a  2-to-l  or  3-to-l  time  advantage  for 
the  minicomputer.  However,  when  more  complex  instructions 
such  as  multiply  and  divide  must  be  execuited,  the  mini¬ 
computer  has  a  considerable  advantage.  For  example,  for  a 
multiply  of  8  bits  by  16  bits  the  Intel  8080  would  take  an 
estimated  230  microseconds  compared  to  just  over  5.5  micro¬ 
seconds  for  the  Inter-Data  50.  Division  times  are  slightly 
longer,  on  the  order  of  270  microseconds  versus  about  10 
microseconds  for  the  Inter-Data  50. 

These  examples  show  that  considerable  care  must  be  exercised 
In  choosing  a  particular  microprocessor  or  microcomputer 
for  the  manpack,  and  that  speeds  must  be  compared  on  the 
basis  of  typical  complex  operations,  rather  than  just  on 
cycle  time.  In  addition  to  computation  speed  one  must 
carefully  examine  the  type  of  software  that  is  available 


3-4 


for  the  computer  and  how  the  programmer  can  use  this  soft¬ 
ware  to  best  advantage.  Some  of  the  considerations  in 
software  design  are  word  size,  both  for  data  and  instruc¬ 
tions,  the  ability  to  address  a  large  number  of  memory 
locations  (which  influences  the  number  of  words  in  the  pro¬ 
gram),  the  number  of  registers  which  are  available  for 
arithmetic  for  indexing  and  for  other  uses  such  as  push 
down  stacks  and  the  number  and  type  of  addressing  modes 
that  can  be  used  such  as  immediate,  relative,  and  relative 
to  program  counter. 

A  second  major  consideration  is  the  hardware  design  inclu¬ 
ding  the  type  of  technology  used  whether  it  be  P-MOS,  N-MOS, 
C-MOS,  or  bi-polar,  which  will  in  turn  influence  power  dissi¬ 
pation,  voltages,  speed,  compatibility  with  other  parts 
which  determines  the  amount  of  interfacing  required,  and 
finally  the  size  of  the  overall  package. 

The  third  major  consideration  is  how  the  hardware  and  soft¬ 
ware  design  influence  the  system.  Here  one  must  look  at  the 
number  of  external  interrupts,  how  many  chips  are  required 
to  implement  the  microprocessor,  whether  or  not  it  is  micro - 
programmable,  whether  it  has  DMA  capability  if  large  amounts 
of  external  data  must  be  transferred  into  memory,  and  the 
amount  of  software  support  available  from  the  manufacturer, 
including  documentation,  application  notes,  design  aids  and 
field  service.  Other  system  considerations  are  the  price, 
the  availability,  and  whether  or  not  the  part  is  second 
sourced  by  other  manufacturers. 


Software  Support 


The  available  software  support  for  microcomputers  is  very 
important  since  having  to  write  the  software  involves  a  con¬ 
siderable  amount  of  money.  Some  of  the  types  of  software  that 
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should  be  available  with  any  microprocessor  are  an  assembler, 
a  monitor  for  debugging  the  microprocessor,  an  editor  for 
changing  programs  and  a  simulator  that  will  allow  the  micro¬ 
processor  program  to  be  executed  on  a  larger  computer. 

The  ability  of  the  operator  to  program  the  computer  using 
higher  level  languages  such  as  FORTRAN  are  not  as  important 
in  the  manpack  application  as  they  might  appear.  One  reason 
is  that  the  manpack  processor  is  a  high  volume  design,  thus 
any  reduction  in  memory  size  resulting  from  the  generally  more 
efficient  code  produced  by  a  good  programmer  writing  directly 
in  assembly  language  would  probably  result  in  a  cost  savings 
(as  well  as  an  increase  in  speed).  Another  point  to  be  con¬ 
sidered  is  that  a  high  level  language  compiler  can  always  be 
written,  if  desired,  and  through  the  use  of  specialized  micro¬ 
programs  built  into  the  machine  can  embody  features  specifi¬ 
cally  optimized  for  the  manpack  application  such  as  trigonometric 
functions,  matrix  manipulation,  interrupt  handling,  floating 
point  and  data  formatting.  Matrix  manipulations,  in  particular, 
would  be  considerably  faster  if  specialized  instructions  and 
microprogramming  is  available.  A  final  point  against  providing 
a  resident  FORTRAN-like  compiler  is  that  a  high  level  language 
would  be  of  little  use  to  the  manpack  operator  since  he  would 
not  be  doing  any  programming.  If  a  compiler  is  used  at  all  it 
would  be  a  cross-compiler,  i.e.,  a  compiler  running  on  a  larger 
machine  which  produces  machine  code  for  the  manpack  processor, 
probably  in  the  form  of  a  pre-programmed  read -only -memory  (ROM) 
which  would  be  inserted  into  the  processor  as  a  maintenance 
operation.  Also,  for  the  operator  to  be  able  to  program  the 
machine  equipment  for  reading  either  paper  tape  or  magnetic  tape 
would  have  to  be  provided  for  program  entry  and  these  are 
definitely  unsuitable  for  the  environmental  conditions  under 
which  the  manpack  will  be  used. 
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It  is  likely  that  one  or  more  standard  microprocessors  will 
be  chosen  for  military  applications  in  order  to  reduce  inven¬ 
tory  requirements  and  to  have  some  commonality  in  software. 

At  the  present  time  there  is  no  clear  picture  of  what  will 
actually  be  available  by  1980  for  military  application  since, 
most  of  the  devices  are  intended  for  civilian  applications. 
However,  there  does  not  seem  to  be  any  doubt  that  sufficient 
technology  will  exist  by  1980,  given  that  currently  available 
microprocessors  and  semi-conductor  memories  could  in  fact  be 
used  with  some  compromises.  Rather,  it  is  more  a  question 
of  guiding  development  to  produce  a  suitable  product  and  hope¬ 
fully  one  with  a  wide  range  of  application. 

3.5  Available  Microcomputers 

3.5.1  Intel  8080 


The  most  prominent  of  the  8-bit  microcomputers  is  the  Intel  8080 
which  in  many  ways  is  the  current  standard  for  8-bit  micropro¬ 
cessors.  The  present  cost  is  $360  in  single  units  and  is  second 
sourced  by  Advanced  Micro  Devices.  The  chip  is  an  N-MOS  dynamic 
technology  and  uses  +  5  volt  and  12  volt  power  supplies.  It 
provides  an  extensive  software  instruction  set  including  data 
manipulation  instructions  for  arithmetic  and  logic,  and  has 
both  BCD  arithmetic  and  double  precision  instructions  (to  string 
two  8-bit  data  bits  together  as  a  16-bit  word).  Data  movement 
instructions  use  three  pairs  of  general  purpose  registers  as 
memory  address  pointers  to  address  both  low  and  high  order 
bits  of  16-bit  memory  addresses  and  can  do  multiple  indexing 
at  the  expense  of  additional  steps.  It  also  has  a  stack  pointer 
which  can  be  used  to  create  push-down  stacks  in  external  memory 
which  would  allow  unlimited  sub-routine  nesting.  It  also  has  a 
multiple  interrupt  capability  and  can  have  direct  memcry  access. 
Intel  provides  a  high  level  language  which  they  call  PL/M,  which 
is  based  on  the  Pl/1  language. 


3-7 


3.5.2 


Motorola  MC  6800 


An  8-bit  microprocessor  made  by  Motorola  (the  MC  6800)  is  also 
a  potential  candidate.  The  present  cost  is  $360  for  one  unit 
but  it  is  estimated  the  price  will  drop  to  $30  in  large  quan¬ 
tities  by  the  end  of  1975.  American  Micro  Systems  also  second 
sources  this  unit.  The  architecture  of  the  machine  is  very 
similar  to  that  of  the  Digital  Equipment  Corporation's  PDP-11, 
except  that  it  has  an  8-bit  word  length.  Execution  times  for 
typical  instructions  are  2-5  microseconds.  The  chip  is  built 
using  N-M0S  silicon  gate  technology  and  requires  only  a  single 
+  5  volt  supply.  One  of  the  interesting  features  of  the 
Motorola  device  is  that  a  large  number  of  external  chips  are 
available  which  interface  directly  to  the  microprocessor 
greatly  simplifying  the  system  design. 

As  mentioned  previously,  the  software  is  very  close  to  that  of 
the  PDP-11  and  provides  data  manipulation  instructions  including 
both  arithmetic  and  logic  as  well  as  instructions  to  take  advan¬ 
tage  of  the  presence  of  two  accumulators  on  the  chip.  Data 
movement  instructions  include  instructions  for  reaching  the 
first  256  memory  locations  or  base  page  with  short  instructions 
and  it  can  process  lists  efficiently  using  the  index  register 
and  relative  addressing.  The  program  manipulation  instructions 
have  most  of  the  features  of  the  PDP-11  software  including 
branches  and  conditional  branches,  as  well  as  unlimited  sub¬ 
routine  nesting  by  using  e  stack  pointer  to  address  a  push  down 
stack  in  external  memory.  One  feature  it  does  not  have  is 
vectored  interrupt  but  it  can  achieve  this  function  at  some 
additional  expense  in  software.  Instructions  for  storing  the 
status  register  are  also  provided  which  is  useful  when  processing 
interrupts.  The  software  support  that  is  currently  available 
is  a  cross-assembler ,  an  interactive  simulator  and  a  monitor 
debug.  In  addition,  high  level  language  based  or  similar  to 
National  Semiconductor's  PL/M  is  currently  being  written.  65K 
8-bit  bytes  of  memory  can  be  addressed. 
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3.5.3 


RCA  COSMAC 


Another  8-bit  microprocessor  with  characteristics  that  would 
lend  itself  to  manpack  use  is  made  by  RCA.  Unlike  the  pre¬ 
viously  mentioned  devices  this  chip  is  based  on  complementary 
MOS  (C-MOS)  which  gives  it  considerably  lower  power  consumption, 
improved  noise  immunity  and  the  ability  to  operate  over  the 
military  temperature  range.  This  chip  will  be  second  sourced 
by  American  Memory  Systems  and  will  be  initially  supplied  as 
two  chips  with  a  single  chip  version  expected  in  late  1975. 

The  total  power  consumption  is  less  than  100  milliwatts.  Data 
manipulations  instructions  include  add,  subtract,  logical, 
reverse-subtract  (which  subtracts  the  contents  of  the  accumu¬ 
lator  from  the  contents  of  a  memory  location),  1-bit  right 
shifts  and  the  ability  to  do  "immediate  instructions."  Data 
movement  instructions  allow  the  relatively  simple  addressing 

of  up  to  65K  of  memory  by  using  a  4-bit  register  to  select 
pointers  in  a  16  x  16  array  which  then  defines  the  memory  loca¬ 
tion  to  be  operated  on.  The  16  x  16  array  also  gives  the  micro¬ 
processor  the  ability  to  jump  to  sub-routines  by  having  a  4-bit 
register  address  different  registers  in  the  16  x  16  array  as  if 
they  were  the  program  counter.  Return  from  sub -routines  is 
done  in  a  similar  manner  by  readdressing  the  original  general 
purpose  register  as  the  program  counter.  Additional  stacks 
can  be  created  in  external  random  access  memory  by  using  the 
pointers  to  index  these  stacks.  Software  support  available  at 
the  present  time  is  in  a  cross-assembler  written  in  Fortran  IV 
which  is  available  on  time-sharing  networks. 
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3.5.4 


General  Instrument  CP-1600 


One  of  the  disadvantages  of  8-bit  architecture  is  the  additional 
time  required  to  execute  instructions  that  require  the  manipula¬ 
tion  of  more  than  8  bits  of  data.  For  this  reason  machines 
with. 12  and  16-bit  capacities  appear  to  be  preferred.  One 
machine  which  is  an  interesting  hybrid  between  an  8 -bit  and  a 
16-bit  processor  is  the  CP-1600  made  by  General  Instrument 
Corporation.  Although  the  arithmetic  logic  unit  is  an  8-bit 
device  the  machine  externally  appears  to  be  a  16-bit  machine. 

The  architecture  of  this  machine  is  similar  (like  many  of  the 
others)  to  a  PDP-11.  The  typical  instruction  speed  execution 
time  is  1.6  to  4.8  microseconds.  Data  manipulation  instructions 
that  are  available  include  add,  subtract  and  logicals,  a  double 
precision  add,  subtract  and  logicals  which  could  be  very  useful 
in  the  high  precision  navigation  routines,  as  well  as  arithmetic 
and  logical  shifts  of  1  and  2  bits.  It  is  also  similar  to  the 
PDP-11  in  that  eight  of  the  additional  general  purpose  registers 
on  the  CPU  chip  can  be  used  as  additional  accumulators.  It  can 
also  address  1024  words  in  the  base  page  and  like  the  PDP-11, 
treats  all  peripheral  devices  as  if  they  were  memory.  The  pro¬ 
gram  manipulation  instructions  available  are  conditional  branches 
and  the  ability  to  address  memory  locations  relative  to  the 
program  counter  up  to  +1024  words  is  provided.  The  machine  also 
has  priority  interrupt  with  self-identifying  vectors  which  will 
greatly  speed  up  the  response  time  to  interrupts.  One  slight 
disadvantage  is  that  no  index  register  is  available  but  the 
same  function  can  be  achieved  using  one  of  the  general  purpose 
registers  as  the  index  register.  Available  software  includes 
a  cross-assembler  and  debug  monitor,  simulator,  diagnostics, 
as  well  as  a  sub-routine  library. 
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IMP -8 /IMP-1 6 


The  IMP- 8  and  IMP- 1.6  are  8  and  16  bit  multichip  microprocessors 
manufactured  by  National  Semiconductor,  These  microprocessors 
ar.e  based  on  a  4-bit  CPU  slice  called  a  RALU  which  has  7 
registers  (4  of  which  are  directly  accessible  to  the  programmer), 
a  16  word  push-down  stack  and  a  four-bit  arithmetic  logic  unit 
(ALU).  Instructions  are  decoded  by  the  Control  Readonly 
Memory  (CROM)  and  used  to  control  the  RALU  up  to  100  23  bit 
micro-instructions  are  stored  in  each  CROM.  Additional  CROM 
can  be  paralleled  to  add  user-supplied  instructions  or  instruc¬ 
tion  set  extensions. 

One  CROM  presently  available  extends  I/O  instructions  to  permit 
block  transfer  of  data,  memory  search  and  stack-store/restore. 
Another  available  CROM  includes  extensions  for  multiply-divide , 
double  precision  add  suucrac  t,  bit  testing  and  byte-handling. 
Typical  instruction  times  are  4.6  nsec  for  a  register-to- 
register  add,  7.7  nsec  for  a  memory-to-register  add  and  10.5 
nsec  for  a  register  input-output  instruction. 

The  typical  instruction  is  executed  in  10  nsec  with  a  2  nsec 
microcycle  time  (i.e.,  4  to  5  microcycles/instruction). 

Complex  operations  involving  many  instructions  take  considerably 
longer.  For  example  software  multiplication  of  a  16  bit  number 
in  one  accumulator  by  a  16  bit  number  in  another  accumulator 
to  produce  a  32  bit  product  requires  almost  one  millisecond. 
However,  in  comparing  these  execution  times  with  those  for 
eight  bit  machines  it  must  be  remembered  that  the  precision 
here  is  twice  that  of  the  eight  bit  machines. 

A  number  of  system  interface  chips  are  being  developed  to 
reduce  the  total  system  package  count.  In  addition^  a  stand¬ 
alone  development  system  with  CPU,  RAM,  ROM  and  I/O  interface 
card  is  available.  National  also  plans  to  introduce  a  high 
level  language  based  on  PL/1  called  SM/PL. 
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3.5.6 


PACE 


The  PACE  microprocessor  by  National  Semiconductor  is  a  16 
bit  single  chip  version  of  their  earlier  multichip  IMP-16. 

It  has  45  instructions,  four  general  purpose  accumulators  that 
are  accessible  to  the  programmer,  a  10  word  push-down  16  bit 
stack  (LIFO)  and  six  vectored  priority-interrupt  levels.  The 
push-down  stack  is  particularly  useful  for  nesting  the  return 
addresses  during  subroutine  branching.  The  general  purpose 
accumulators  give  a  great  deal  of  flexibility  in  addressing 
modes.  PACE  can  also  handle  8  bit  data  efficiently.  Both 
direct  and  indirect  memory  addressing  are  provided.  Direct 
memory  addressing  permits  the  programmer  to  specify  a  memory 
location  that  is  either  on  the  base  page  (256  words)  or  +  127 
locations  relative  to  either  the  program  counter  or  index 
register.  The  relative  addressing  helps  to  prevent  the  rather 
small  base  page  from  being  used  entirely  just  by  linkages 
between  pages.  Indirect  addressing,  which  uses  the  contents 
of  a  particular  memory  location  as  the  address,  can  be  done 
relative  to  two  of  the  four  accumulators  used  as  index  registers, 
relative  to  the  program  counter,  or  using  the  base  page.  Up 
to  65K  words  of  memory  can  be  addressed  in  this  manner.  Since 
I/O  devices  are  treated  as  memory  locations  the  available 
addresses  may  be  divided  between  RAM,  ROM  and  I/O  devices  as 
desired. 
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3.6  Microcomputer  System  Manufacturers 

A  number  of  companies  offer  microprocessor-based  systems  built 
aro.und  either  their  own  or  someone  else's  chips.  These 
systems  are  useful  in  the  developmental  and  prototype  stage; 
moreover,  some  of  these  systems  include  proprietary  software 
not  available  from  the  chip  manufacturer.  Table  3-1  (Falk, 
1974)  lists  manufacturers  and  characteristics  of  a  number  of 
available  microcomputer  development  systems. 

Table  3-2  shows  the  relative  speed  and  memory  requirements  of 
some  of  the  more  widely  available  processors  when  running 
typcial  small  tasks.  These  numbers  are  plotted  in  Figure  8. 
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4.0 


COMPUTER  FUNCTIONS 


4.1  Introduction 


This  chapter  describes  the  major  functions  which  the  computer 
must  perform  including  both  computation  and  receiver  control. 

The  topic  of  recursive  filters  for  position  estimation  and 
its  implications  for  both  computer  selection  and  system  per¬ 
formance  is  discussed.  Of  particular  interest  are  alternate 
forms  which  are  less  sensitive  to  computer  word  lengths  and 
thus  are  more  suitable  for  implementation  in  small  computers. 
Other  topics  include  estimates  of  computation  time  as  a  func¬ 
tion  of  the  number  of  filter  states,  causes  of  filter  diver¬ 
gence  and  initialization  procedures. 

The  above  topics  are  motivated  by  the  need  to  use  small  micro¬ 
processors  in  the  manpack  in  order  to  meet  size,  weight,  power 
consumption  and  cost  goals. 

4.2  Receiver -Computer  Interactions 

In  addition  to  its  computational  functions  the  manpack  com¬ 
puter  will  also  control  and  sequence  receiver  operations  asso¬ 
ciated  with  both  initial  satellite  acquisition  and  housekeeping. 
To  illustrate  some  of  the  tasks  which  the  computer  must  perform 
the  sequence  of  operations  required  to  utilize  the  P-code  are 
listed  below: 

1.  Select  constellation  based  on  approximate  location 
and  time  using  Keplerian  Orbit  model  and  GDOP 
calculation. 

2.  Compute  expected  doppler  and  set  doppler  offset 
into  receiver. 
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3. 


Select  proper  code,  load  into  receiver  code 
generator  and  step  generator  in  one-half  chip 
increments  until  lock  is  achieved  or  some  predeter¬ 
mined  length  of  time  passes  without  lock  being 
obtained.  In  the  latter  case  return  to  the  error¬ 
handling  routine. 

4.  When  lock  is  achieved  locate  the  sync  and  Hand- 
Over-Word  (HOW)  in  the  data. 

5.  Using  the  sync  and  HOW,  acquire  the  P-code  and 
measure  the  pseudorange. 

6.  Demodulate  the  P-code  data  for  fine  orbit  position 
determination. 

Steps  2-6  are  performed  in  parallel  for  4-channel 
receiver  and  sequentially  for  sequential  receiver. 

7.  Using  recursive  filtering  compute  the  receiver 
position  and  clock-bias  based  on  the  fine  orbit 
model  with  corrections  for  ionospheric  and  tropo¬ 
spheric  delay  (using  L1/L2  corrections  or  models). 

8.  Periodically  check  constellation  for  GDOP  and 
elevation  angle.  If  better  constellation  becomes 
available  change  to  new  constellation.  If  con¬ 
stellation  deteriorates  badly  and  no  new  constella¬ 
tion  is  available  go  to  error -hand ling  routine. 

9.  Transform  receiver  position  into  UTM  coordinates 
as  required. 

\ 

Although  some  of  the  operations  listed  above  could  be  done 
either  in  software  or  hardware  they  have  been  partitioned 
such  that  the  computer  performs  those  operations  which  either 
require  considerable  computation  or  would  require  special 
purpose  hardware  with  a  low  utilization.  These  tasks  include 
filtering,  constellation  selection,  doppler  prediction,  pseudo¬ 
range  and  pseudo-doppler  calculation  and  data  demodulation. 


Hardware  is  most  useful  for  implementing  high  speed  or  repeti¬ 
tive  tasks  such  as  code  generation,  acquisition  and  tracking. 
System  time  also  should  be  maintained  with  hardware.  Either 
software  or  hardware  could  be  used  to  locate  the  sync  and 
HOW.  In  any  case  the  computer  must  have  data  such  as  the 
coarse  orbital  model  parameters,  code  words  for  each  satellite 
and  ionospheric  correction  model  parameters  stored  in  read- 
only-memory  (ROM)  to  prevent  loss  of  this  data  when  the  equip¬ 
ment  is  turned  off. 

Because  of  the  number  of  computer -receiver  interactions  during 
the  initial  acquisiton  process  the  manpack  processor  must  have 
an  efficient  interrupt-handling  scheme.  Polling,  where  the 
computer  continually  loops  through  a  segment  of  code  to  check 
the  status  of  external  hardware,  would  waste  computer  capa¬ 
bility  in  this  application. 

One  real  advantage  of  incorporating  as  much  of  the  receiver 
control  as  feasible  in  software  is  that  it  becomes  possible 
to  modify  receiver  operation  without  extensive  hardware  modi¬ 
fications,  which  may  be  particularly  useful  in  early-generation 
receivers.  In  addition,  the  smaller  amount  of  hardware  reduces 
initial  cost  and  power  consumption,  increases  reliability  and 
simplifies  maintenance.  From  a  systems  point -of-view  these 
benefits  could  be  well  worth  any  reduction  in  speed. 

4.3  User  Position  Determination 

In  order  for  the  user  to  compute  position,  the  receiver  must 
first  accurately  measure  the  corrected  pseduorange  to  each 
spacecraft.  The  pseudorange  measurement  in  conjunction  with 
the  accurately  known  spacecraft  orbital  position  allows  the 
computation  of  position;  however,  system  time  must  also  be 
accurately  known  in  order  to  do  the  precise  orbit  determination. 
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This  time  information  (user  clock  bias)  can  be  obtained  as 
one  of  the  solve-  -variables  in  a  recursive  estimation 
procedure.  The  accurate  system  time  thereby  derived  is  used 
in  a  precision  orbit  position  model  along  with  the  orbit  model 
parameters  transmitted  as  data. 


4.3.1 


Recursive  Filterim 


Recursive  filtering  is  used  extensively  in  navigation  appli¬ 
cations  since  it  allows  real-time  position  updating  by  incor¬ 
porating  new  measurement  data  while  at  the  same  time  reducing 
the  error  by  filtering.  The  best  known  of  these  techniques  is 
Kalman  filtering  and  its  computational  variations  such  as  the 
Joseph  measurement  incorporation  algorithm  and  the  Potter 
square  root  filter.  These  filters  are  of  great  interest  be¬ 
cause  they  are  capable  of  providing  the  minimum  mean  square 
linear  estimate  of  a  set  of  parameters  corrupted  by  noise. 

The  equations  for  these  filters  are  presented  in  Appendix  C, 


Although  the  Kalman  formulation  is  efficient  computationally, 
if  the  problem  is  ill-conditioned  the  measurement  incorporation 
algorithm  of  the  standard  Kalman  filter  is  very  sensitive  to 
computer  word  length  round-off.  The  Potter  square-root  design 
has  approximately  twice  the  precision  of  the  Kalman  form  for 
ill-conditioned  problems  and  is  thus  more  stable  (Kaminski, 
1971);  however,  this  increase  in  stability  is  obtained  at  the 
expense  of  significantly  greater  computation  time.  Carlson 
(Carlson  1973)  has  developed  a  modified  form  of  the  square 
root  filter  which  increases  computation  speed  to  approximately 
that  of  the  standard  Kalman  filter  in  low-order  designs  while 
retaining  the  numerical  stability  of  the  square  root  imple¬ 
mentation.  Carlson's  paper  contains  tables  which  compare  the 
number  and  typa  of  arithmetic  operations  needed  for  measure¬ 
ment  incorporation  and  time  update  steps  in  these  different 
filter  implementations.  These  tables  are  reproduced  as  Tables 
4-2  and  4-3. 
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Table  4,2  Number  of  arithmetic  operations  for  one 
measurement  incorporation  step 
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Table  4.3  Number  of  arithmetic  operation  for 
one  time  update  step 
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In  these  tables  m  is  the  number  of  simultaneous  measurements 
to  be  incorporated  in  each  measurement  step,  n  is  the  number 
of  filter  states,  n^  is  the  number  of  independent  process 
noises  being  considered  and  RSS  and  Householder  refer  to 
alternative  techniques  for  triangularization  of  the  updated 
covariance  square-root  matrix.  The  Householder  triangulariza¬ 
tion  requires  about  one-half  the  numerical  precision  of  the  RSS 
method  but  needs  appreciably  more  computation  time. 


4.3.2 


Effect  of  the  Number  of  Filter  States  on  Pro¬ 
cessor  Design  and  Filter  Performance  ~ * 


Tables  4-2  and  4-3  show  that  the  number  of  computations  rises 
rapidly  as  the  number  of  states  in  the  filter  increases.  For 
a  given  processor  this  increase  in  the  number  of  required 
computations  causes  the  minimum  time  between  updates  to  increase, 
with  the  possible  result  being  poorer  overall  performance  since 
measurements  are  not  being  incorporated  at  as  high  a  rate  as 
perhaps  is  possible,  causing  a  loss  in  information.  For  the 
moving  user,  in  particular,  excessively  long  intervals  between 
updates  will  greatly  increase  the  error.  Conversely,  if  a 
sufficient  number  of  states  is  not  included  in  the  filter 
model  then  the  filter  performance  will  also  be  poor  and  can 
in  fact  diverge.  Clearly,  then,  there  is  an  optimum  number  of 
filter  states  based  on  the  user  dynamics,  noise  sources, 
available  computation  speed  and  accuracy,  number  of  measure¬ 
ments  and  filter  algorithm  used.  Ultimately,  candidate  filters 
must  be  tested  in  actual  equipment  or  with  an  accurate  simula¬ 
tion  which  Includes  computer -dependent  effects  and  all  of  the 
sources  of  error.  A  covariance  analysis  is  not  sufficient 
since  it  does  not  include  computational  effects.  One  possible 
result  of  failure  to  test  a  filter  design  under  conditions 
that  adequately  simulate  actual  use  is  a  design  that  performs 
substantially  poorer  than  predicted. 
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Under  some  conditions  the  filter  performance  error  can  become 
so  large  that  the  estimates  are  useless.  When  this  happens 
the  filter  is  said  to  have  diverged.  This  behavior  is  typically 
caused  by  a  poor  model  of  the  error  sources  which  results  in 
a  covariance  matrix  which  is  too  small  (and  thus  a  higher  gain 
filter)  and  does  not  properly  weight  the  measurements.  Diver¬ 
gence  may  result  from  one  or  more  of  the  following: 

1.  Truncation  errors  caused  by  finite  word  length. 

2.  Failure  to  model  the  noise  sources  properly,  i.e 
poor  estimates  of  the  state  error  covariance  matrix, 
the  process  noise  covatiance  matrix,  or  the  measure¬ 
ment  error  covariance  matrix.  Failure  to  include  a 
significant  noise  source  can  be  disasterous. 

3.  Failure  to  model  the  process  itself  properly. 

4.  Failure  to  model  the  measurement  process  correctly. 

Since  one  of  the  possible  causes  of  filter  divergence  is  the 
failure  to  accurately  model  or  include  noise  sources  it  is 
important  that  the  designer  identify  all  of  the  significant 
contributors  along  with  estimates  of  their  magnitudes.  These 
sources  include  the  following: 

1.  Lack  of  knowledge  of  true  system  time  because  of 
drifts  in  the  phase  and  frequency  of  the  satellite 
oscillator  which  are  not  accurately  modeled.  Both 
random  and  bias  errors  are  present. 

2.  Inability  to  measure  the  pseudorange  accurately  as 

a  result  of  propagation  delays  that  are  not  accounted 
for  by  the  model.  These  effects  will  be  increasingly 
important  at  low  elevation  angles,  and,  unless  L1/L2 
measurements  are  made,  will  contribute  both  random 
and  bias  errors. 
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3.  Errors  in  measurement  of  code  phase  due  to  the  pre¬ 
sence  of  thermal  and  jamming  noise  in  the  receiver 
and  the  effect  of  code  resolution.  These  errors 
are  zero  mean. 

4.  Errors  in  estimating  doppler  (if  used  as  part  of 
the  filter  states)  depend  on  the  measurement  tech¬ 
nique  but  in  general  have  similar  sources  as  3. 

5.  Inability  to  accurately  know  the  position  of  the 
satellite  at  a  given  time.  Obviously  the  satellite 
position  (ephemeris)  must  be  known  accurately  to 
compute  position  and  user  time.  Error  sources  are 
the  satellite  position  model  stored  on  the  satellite 
and  satellite  clock  errors. 

6.  Errors  in  pseudorange  and  user  time  bias  calcula¬ 
tions  resulting  from  drift  and  noise  in  the  user 
oscillator. 

Most  techniques  that  have  been  proposed  for  avoiding  filter 
divergence  are  somewhat  empirical  and  must  be  adjusted  accord¬ 
ing  to  the  application.  One  of  these,  "over  weighting'* 
(Schmidt,  1967)  allows  a  weighting  factor  that  weight c  the 
measurements  more  heavily  than  the  standard  filter,  thereby 
preventing  computationally  small  error  covariances. 

Another  weighting  technique  considered  by  both  Tarn  and 
Zaborsky  (Tarn  and  Zaborsky  1970)  and  Sorenson  (Sorenson 
etal  1970)  is  age  weighting  which  exponentially  weights  older 
data  such  that  new  data  is  emphasized,  thereby  overweighting 
newer  data.  It  has  the  advantage  of  being  simple  to  imple¬ 
ment  with  either  covariance  or  information  filters  by  simply 
replacing  the  state  transition  matrix  ^  by  where 

^k*  "  •  exp  03Ak), 
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P  >  0,  and  A^  is  the  time  between  stage  k  and  k+1.  Thus  A^ 

is  the  "time -constant"  over  which  most  of  the  filtering  is 
done. 

In  general  these  techniques  should  only  be  applied  if  more 
straight-forward  techniques  such  as  square-root  filters  fail 
since  they  compromise  the  filter  performance  if  used  unnecessary- 
ily. 


4.3.3  Recursive  Filter  Computation  Estimates 

The  minimum  number  of  states  is  application  dependent;  however 
to  provide  a  basis  for  estimating  computational  requirements 
three  possible  classes  of  users  are  proposed: 

1.  Stationary  User  -  filter  states  are  the  user 
oscillator  phase  and  frequency  offset  plus  user 
position  (5  states).  Pseudorange  (PR)  and 
pseudo-doppler  (PD)  are  measured. 

2.  Moving  User  I  -  filter  states  are  the  user 
oscillator  phase  and  frequency  offset  plus  both 
user  position  and  velocity  (8  states).  Both  PR 
and  PD  measured. 

3.  Moving  User  II  -  the  filter  states  include  those 
of  Moving  User  I  with  the  addition  of  user 
acceleration  (11  states).  Both  PR  and  PD 
measured. 

The  number  of  additions,  multiplications,  divisions  and  square 
roots  required  to  implement  filters  for  each  of  the  above 
three  classes  were  calculated  using  the  formulas  in  Tables 
4-2  and  4-3.  The  amount  of  computation  time  required  for  the 
three  classes  was  estimated  by  using  published  average  exe¬ 
cution  times  for  the  Hewlett  Packard  Model  2116  computer 
library  routines  shown  in  Table  4-4. 
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Table  4-4 


Execution  Time  for  Hewlett  Packard  2116 
Computer  Library  Routines  (1.6  microsecond  cycle  time) 


Execution  Time  (Milliseconds) 


Routine  Name 

Minimum 

Average 

Maximum 

Floating  Point  Addition 

0.3 

0.5 

0.9 

Floating  Point  Multiply 

0.64 

0.7 

0.75 

Floating  Point  Divide 

1.2 

1.3 

1.5 

Floating  Point  Square  Root 

3.94 

A. 8 

6.45 

These  execution  times  for  the  Hewlett  Packard  2116  are  typical 
of  16  bit  minicomputers  and  some  of  the  faster  microcomputers. 
These  estimates  are  only  approximate  since  there  will  be  vari¬ 
ations  resulting  from  cycle  time  and  architectural  differences. 
These  times  also  assume  that  the  filter  is  implemented  with 
software  floating  point  computations  and  that  the  resulting 
accuracy  of  approximately  24  bits  including  sign  is  sufficient 
(floating  point  numbers  are  represented  by  two  16-bit  computer 
words  with  the  exponent  and  its  sign  occupying  eight  bits  and 
the  fraction  and  its  sign  occupying  24  bits).  Inclusion  of  a 
hardware  floating  point  processor  would  greatly  increase  the 
filter  speed. 

The  computation  requirements  and  estimated  time  for  one 
measurement  incorporation  and  one  time  update  for  the  three 
classes  of  users  are  shown  in  Tables  4-5  through  4-7.  The 
numbers  in  parentheses  are  the  times  required  for  each  type  of 
arithmetic  operation  based  on  Table  4-4.  These  times  do  not 
include  the  operations  required  for  computing  the  measurement 
noise  covariance  matrix,  the  measurement  gradient  matrix  or 
the  measurement  residuals.  It  is  also  assumed  (Carlson  1973) 
that  the  measurement  noise  processes  are  uncorrelated  and  that 
measurements  are  incorporated  one  at  a  time,  except  for  the 
Joseph  square  root  filter.  These  estimated  times  do  not  in¬ 
clude  overhead  time  for  computing  subscripts  for  matrix 
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Total  Time  (Sec)  0.887  1.137  1.156  1.438  0.847 
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operations  (which  only  require  the  much  faster  single  pre¬ 
cision  integer  arithmetic)  nor  do  they  account  for  possible 
inefficient  code  generation  by  a  high  level  compiler  relative 
to  assembly  language  coding. 

Examination  of  the  estimated  times  for  each  of  the  possible 
user  classes  shows  that  the  Carlson  square-root  filter  and 
the  Kalman  filter  are  clearly  superior.  However  the  better 
stability  of  the  Carlson  square-root  filter  in  conjunction 
with  Householder  triangularization  make  it  the  recommended 
technique. 

The  estimated  time  for  one  measurement  incorporation  plus  one 
filter  update  for  each  user  class  with  the  Carlson  square  root 
filter  and  Householder  triangularization  is: 

a.  Stationary  User  -  0.761  seconds 

b.  Moving  User  I  -  1.310  seconds 

c.  Moving  User  II  -  3.529  seconds 


For  a  complete  constellation  of  four  satellites  and  a  four 
channel  receiver  these  times  would  be  increased  by  a  factor 
of  four. 

For  the  single  channel  receiver  these  times  must  include  any 
additional  time  for  reacquiring  the  next  satellite  in  the 
constellation.  Whether  this  process  actually  results  in 
increased  time  depends  on  the  relative  duration  of  the 
measurement  interval  to  the  filter  measurement  incorporation , 
reacquisition,  and  filter  update  time.  It  is  possible  to 
allow  the  computer  to  be  processing  the  most  recent  measure¬ 
ments  while  the  receiver  is  reacquiring  the  next  satellite  and 
making  the  next  measurement.  Whichever  process  is  the  longest 
will  set  the  minimum  filter  cycle  time. 
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Normally,  the  receiver  is  capable  of  making  measurements  at  a 
rate  in  excess  of  the  ability  of  the  filter  to  incorporate 
them,  particularly  in  the  case  of  a  four-channel  receiver. 

This  mismatch  means  that  the  filter  cannot  make  use  of  all 
of  the  available  information,  leading  to  poor  performance. 
Joglekar  (Joglekar,  1973)  suggested  a  technique  for  incorporating 
these  additional  measurements  which  consists  of  simply 
averaging  several  successive  measurements  and  using  the 
average  in  the  measurement  incorporation  step.  Obviously,  the 
number  of  measurements  that  can  be  handled  in  this  way  is 
limited  by  the  user  dynamics  although  typically  2  to  4 
measurements  can  be  averaged.  It  appears  that,  at  least  for 
the  stationary  user,  measurement  averaging  can  be  a  useful 
technique  for  reducing  cede -phase  measurement  errors  resulting 
from  thermal  noise,  jamming  and  oscillator  phase  noise  since 
the  expected  RMS  error  will  decrease  as  (N)^  where  N  is  the 
number  of  measurements  averaged. 

In  practice  there  may  be  25  to  50  cycles  before  the  filter 
estimate  has  sufficiently  small  error  to  be  useful.  The  time 
required  for  these  cycles  must  added  to  the  acquisition 
time  to  estimate  the  length  of  time  before  a  useful  "fix"  is 
obtained.  Unless  computation  speed  can  be  increased  by  tech¬ 
niques  such  as  a  hardware  floating-point  processor  added  to 
the  basic  microprocessor  or  more  efficient  instruction  execu¬ 
tion,  i.e.,  by  microprogramming,  it  would  appear  that  execu¬ 
tion  times  are  marginal  except  for  the  stationary  user.  Micro¬ 
programming  is  preferred  over  a  hardware  floating  point  pro¬ 
cessor  because  of  the  increased  power  and  cost  that  it  would 
entail.  However,  if  even  greater  speed  is  required  the  hard¬ 
ware  processor  is  the  next  best  choice. 

4.3.4  Storage  Requirements  for  Filter  Variables 

Carlson  (Carlson, 1973)  gives  the  program  variable  storage 

requirements  in  Table  A-8.  The  estimated  number  of  storage 

\ 
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locations  for  each  of  the  three  classes  of  user  based  on 
Table  4-8  is  given  in  Table  4-9.  Since  these  are  floating 
point  numbers  occupying  32  bits  each  the  number  of  locations 
identified  in  Table  4-8  should  be  increased  by  a  factor  of  2 
for  a  16  bit  memory  and  by  a  factor  of  4  for  an  8  bit  memory. 


Table  4-8  Program  Variable  Storage  Requirements 


Filter  formulation 

Variable  storage  words 

New  square  root  (RSS) 

n^+n^+2n2+3n 

New  square  root  (Housr.) 

«S23fi  —  ! 

Potter  square  root  (RSS) 

n^ -hi  +2  n2 +3  n 

Potter  square  root  (Housr.) 

n2+p2+2n2+3n 

Kalman  conventional 

n^-hi^^+2n2+3n 

Joseph  conventional 

2n^+2mn+m^+2n 

Table  4-9  Filter  Variable  Storage  Requirements 


Moving 
User  II 

Moving 
User  I 

Stationary 

User 

Carlson  SR  with  RSS 

275 

152 

65 

Carlson  SR  with 
Householder 

275 

152 

65 

Potter  SR  with  RSS 

275 

152 

65 

Potter  SR  with 
Householder 

275 

152 

65 

Kalman 

275 

152 

65 

Joseph 

312 

180 

84 
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As  can  be  seen  from  this  table  the  required  amount  of  storage 
is  the  same  for  all  of  the  various  formulations  with  the 
exception  of  the  Joseph,  which  requires  about  30%  more  storage 
for  the  stationary  user  and  about  13%  for  the  moving  user  II. 
In  any  case  the  required  number  of  storage  locations  is  not 
particularly  large.  It  should  be  noted,  however,  that  unlike 
the  program  storage  which  will  most  likely  be  incorporated 
in  a  read-only  memory  (ROM)  the  program  variable  storage  must 
be  read-write  random-access  memory  (RAM). 

4.3.6  Explicit  Navigation  Solution  , 

In  contrast  to  recursive  (Kalman)  filtering  of  pseudorange 
measurements,  position  also  may  be  calculated  explicitly. 

One  technique  for  explicit  calculation  (Schmidt,  1972)  is 
relatively  simple.  The  calculation  is  based,  in  three  dimen¬ 
sions,  on  finding  the  focus  of  a  three  dimensional  conic 
defined  by  the  four  satellites  and  the  measured  pseudorange 
differences.  It  can  be  shown  that  the  focus  is  the  user 
location. 

The  solution  proceeds  by  selecting  three  of  the  satellites  and 
finding  the  equation  ol:  a  plane  containing  the  focus.  This 
procedure  is  repeated  for  the  remaining  triads  and  the 
resulting  set  of  linear  equations  solved  for  the  common  inter¬ 
section  of  planes,  which  is  the  focus  and  hence  the  user 
location.  The  equation  for  the  plane  is: 

(XjA23  +  ^31  +  X3a12)X  +  (Yj*23  +  Y2431  +  Y3*12)Y 

+  (Zj*23  +  Z2*31  +  Z3*12)Z 

“  !/*  (*12*23*31  +  al  *23  +  a22  *31  +  a32  *12^ 

which  is  of  the  form  Ax  +  By  +  Cz  *  D,  where  is  the  range 
difference  between  the  ith  and  jth  satellite,  and 
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+  Zj^,  where  (Xj,  Y y  Zj)  are  the  coordinates 
of  the  jth  satellite.  Thus 


A1  B1  C1 

X 

Di 

A2  B2  C2 

Y 

- 

D2 

A3  B3  C3 

Z 

D3 

L  m 

and  the  subscripts  identify  the  triad.  It  should  be  noted 
that  four  triads  can  be  chosen  from  four  satellites.  Since 
only  three  triads  provide  independent  solutions,  the  fourth 
triad  can  either  be  ignored  or  included  via  a  least -mean- 
square  solution.  The  latter  approach  is  preferable  in  that 
some  averaging  of  error  results. 

The  least-mean-square  solution  which  incorporates  all  four 
possible  triads  is: 


X 

1 
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A1A2A3A4 
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B1  Cl 

V 

A1A2A3A4 
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D 

< 
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b3  C3 
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C1C2C3C4 

°3 

\ 

A4 

b4  C4 

°4 

If  more  than  four  satellites  are  visible  or  if  additional 
measurements  to  the  same  four  satellites  are  available 
the  above  solution  can  be  extended  to  include  the  additional 
measurements.  One  interesting  aspect  of  the  least-squares 
solution  is  that  the  matrix  to  be  inverted  is  always  3X3. 

Since  this  position  calculation  does  not  account  for  the 
magnitudes  of  the  various  error  sources  and  weight  them 
appropriately  the  error  will  not  be  as  small  as  for  the 
optimum  linear  estimator  (Kalman  filter).  However,  if  the 
magnitudes  of  the  bias  errors  are  larger  than  the  random 
errors  the  Kalman  filter  will  not  perform  well.  Under  these 
conditions  che  easily  calculated  explicit  solution  considered 
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above  may  not  be  much  worse  in  the  case  of  the  stationary 
user.  If  the  user  is  in  rapid  motion  then  the  explicit 
solutions  give  poor  accuracy. 

Explicit  solutions  may  be  useful  in  Kalman  filter  initializa¬ 
tion  since  they  do  not  require  any  a  priori  knowledge  nor  do 
they  require  knowledge  of  user  time  bias.  These  calculations 
can  be  iterated  with  the  improved  estimates  of  system  time 
used  to  improve  satellite  ephemeris  and  hence  the  accuracy  of 
the  fix. 
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EXAMPLE  PRIME  ITEM  DEVELOPMENT  SPECIFICATION 

for  the 

GLOBAL  POSITIONING  SYSTEM 
MANPACK  NAVIGATION  SET 
Type  Bl 


Stanford  Telecomnunications,  Inc. 

1161  San  Antonio  Road 
Mountain  View,  California  94043 


.  ,  .  - -  - - - 


Page  2  of  15 
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Test  Plan 
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1  Jan  74 

DoD  World  Geodetic  System 

Conference,  1972 

Standards 

Federal 

None 

Military 

MIL-STD-454D 

31  Aug  73 
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for  Electronic  Equipment 

MIL-STD-810B 

15  Nov  67 

Change  4 

21  Sep  70 

Environmental  Test  Methods 

MIL- STD- 1472A 

15  May  70 

Human  Engineering  Design  Criteria  for 
Military  Systems  Equipment  &  Facilities 

MIL- STD- 188- 100 

15  Nov  72 

Common  Long  Haul  and  Tactical  Communi¬ 
cation  System  Technical  Standards 

Drawings 
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NTS  PRN  Navigation  Assembly /User  and 
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3.0  Requirements .  The  manpack  set  shall  be  designed  to 

generate  position  and  time  data  when  operated  in  the  GPS  signal 
environment  specified  in  SS-GPS-101.  The  set  shall  be  designed 
to  satisfy  mandatory  requirements  for  site,  weight,  and  contin¬ 
uous*  operating  time  while  optimizing  accuracy,  acquisition  time, 
and  AJ  performance. 

3.1  Item  Definition.  The  manpack  set  shall  consist  of  a 
single  self-contained  package  which  includes  all  major  components 
shown  in  Figure  I  .  It  shall  be  nossible  for  an  operator  to  use 
the  set  in  all  its  normal  operating  modes  while  it  is  being  carried. 

3.1.1  Item  Diagram.  See  Figure  I 

3.1.2  Interface  Definition.  Figure  II  illustrates  the  required 
interfaces  between  the  manpack  and  external  equipment. 

3. 1.2.1  Manpack/Satellite  Interface.  The  Satellite/user  inter¬ 
face  requirements  of  ICD  MH08- 00002- 400  apply. 

3. 1.2. 2  Manpack/RF  Signal  Interface.  The  manpack  shall  provide 
for  connection  to  a  GFE  RF  Signal  Generator  with  characteristics 
specified  in  (TBS). 

3. 1.2.3  Manpack/External  Time  and  Frequency  Source.  The  manpack 
shall  be  capable  of  operation  with  either  internal  or  external 
time  and  frequency  source  which  provides,  as  a  minimum,  one  pulse 
per  second  (1  pps),  one  kilopulse  per  second  (k  lpps),  and  5  MHz 
to  the  set. 

3. 1.2.4  Manpack/External  Power.  The  manpack  shall  be  designed 
to  connect  to  and  operate  on  external  power  for  the  purpose  of 
re-charging  and/or  conserving  internal  battery  power. 


MANPACK  INTERFACES 
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3, 1.2.5  Manpack/Communications  Interface.  The  manpack  shall 
provide  a  synchronous  2.4  kbps  serial  data  and  clock  output  sig¬ 
nal  with  characteristics  as  specified  in  MIL-STD- 188C .  Digital 


data  shall  include  position,  time,  and  status  in  accordance  with 


the  formats  defined  in  Section  10. 


3. 1.2.6  Manpack/Data  Base  Initialization.  The  manpack  design 
shall  provide  for  loading  initial  data  base  parameters  into 
memory  as  required  to  control  signal  acquisition. 


3.1.3  Major  Components.  Major  omponents  of  the  manpack  set  are 

a.  Antenna 

b.  Receiver 

c.  Processor 

d.  Computer  Programs 

e.  Display 

f.  Battery  Pack. 


3.2  Characteristics 

3.2.1  Performance.  The  performance  requirement  of  this  section 

shall  be  satisfied  within  the  physical  constraints  of  size,  and 

weight,  specified  in  3.2.2 

3.2. 1.1  Operating  Modes.  As  a  minimum,  the  set  shall  operate  in 
the  modes  shown  in  Table  1  when  receiving  LI  C/A  channel  signals. 
Mode  selection  shall  be  by  operator  control. 


3.2. 1.2  Operating  Time.  The  manpack  set  shall  include  an  internal 
rechargeable  power  supply  capable  of  operating  for  a  minimum  of  48 
hours  at  32°F  during  which  an  average  of  4  fixes  per  hour  are  made. 
Continuous  operation  for  at  least  12  hours  at  32°F  shall  be  demon¬ 
strated. 
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MODE 

SOLUTION 

PARAMETERS 

OPERATOR 

SPECIFIED 

PARAMETERS 

MINIMUM 

SIGNALS 

RECEIVED 

I 

3- Position 

1-Time 

None 

4 

II 

2- Position 

1-Time 

Alt 

3 

III 

2- Position 

Alt 

Time 

2 

IV 

1  -  Time 

3- Position 

1 

V 

Velocity 

Bearing 

Alt 

None 

4 

TABLE  1 

OPERATING  MODES 
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3 .2 .1.3  Navigation  Error  Performance,  Errors  contributed  to  the 
navigation  solution  by* the  manpack  for  GDOP*s  of  2  or  less  under 
dynamic  and  signal  levels,  specified  herein,  shall  not  exceed: 


a. 

Horizontal  CEP 

10 

meters 

b. 

Vertical  PE 

5 

meters 

c . 

Time  (IT) 

30 

nanoseconds 

excluding  ionospheric  errors.  The  manpack  shall  be  capable  of  re¬ 
ducing  errors  due  to  ionosphere  to  at  least  50%  of  their  nominal 
values. 


3. 2. 1.4  Response  Time  Performance.  The  manpack  design  shall  be 
capable  of  the  following  mutually  exclusive  response  times: 

a.  Preparation  time  prior  to  mission:  10  minutes,  max. 

b.  Time- to- first- fix:  5  minutes,  max. 

c.  Update  rate:  once  per  minute,  max. 

d.  Response  to  manual  request  for  fix:  1  second,  max. 


3.2. 1.5  Dynamics.  Performance  requirements  of  this  paragraph 
shall  apply  under  any  combination  of  the  following  user  dynamic 


conditions: 

Navigation 

Tracking 

a. 

Velocity 

2  meters /sec 

10  meters /sec 

b. 

Acceleration 

2.5  " 

10  " 

c. 

Jerk 

Not  applicable 

Not  applicable 

3.2. 1.6  Signal  Threshold  Levels.  The  manpack  shall  operate  with 
input  signal  levels  from  -130  dbw  to  threshold.  Threshold  levels 
shall  be  defined  as  carrier- to- noise  density  C/No  required  for  the 
following  functions: 

Function  C/No 


Acquisition  Probability  .99  33  db-Hz 
Nav.  Data  Error  Rate  -10" 5  30  db-Hz 
Carrier/Code  Tracking  f'r  =  0.2  27  db-Hz 
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3.2.2  Physical  Characteristics 

3. 2. 2.1  Size /Volume .  1000  cu  inches,  maximum,  including  battery, 

antenna,  control  and  display. 


3. 2. 2. 2  Weight .  12  pounds,  including  internal  battery,  antenna, 

control  and  display. 


3. 2. 2. 3  Power.  Power  consumption  shall  be  consistent  with 

3. 2. 1.2 


3.2.3  Reliability.  The  manpack  shall  possess  a  Mean-Time- 

Be tween- Failure  (MTBF)  of  2000  hours. 

3.2.4  Maintainability.  The  manpack  shall  be  designed  for  a 
Mean-Time-To-Repair  (MTTR)  of  0-10  minutes  at  the  organizational 
level,  20  minutes  direct  support,  and  up  to  45  minutes  at  the 
general  support  level. 

3.2.5  Environment.  (TBS) 

3.2.5. 1  Temperature .  The  manpack  shall  be  capable  of  perform¬ 

ing  to  specification  within  the  temperature  range  of  32°F  to 
+155°F  with  an  ultimate  low  temperature  requirement  of  -40°F. 

3.2. 5.2  Humidity.  Manpack/vehicular  set  performance  shall  not 
degrade  when  subjected  to  relative  humidity  of  100%  at  all  ambient 
air  temperatures  to  85°  and  relative  humidities  corresponding  to 

a  dew  point  of  85°F  at  all  temperatures  between  85°F  and  120°F, 
and  low  relative  humidity  of  5%  at  120°F. 
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3.2.5.3  Elevation. 

a.  Operating  3,000  meters  min.  (-10,000  feet) 

b.  Non-operating  15,000  meters,  min.  (-50,000  feet) 

3. 2. 5. 4  Rain /Immersion. 


3. 2. 5. 5  Shock,  Bench  Handling 


3. 2. 5. 6  Vibration.  The  equipment  shall  be  capable  of  specified 
performance  when  mounted  on  either  wheeled  or  tracked  vehicles.  To 
meet  this  requirement,  the  equipment  may  be  designed  to  operate 
either  with  or  without  vibration  isolators, 

3.2.5. 7  Orientation.  Excepting  antenna  effects,  performance  of 
the  manpack  shall  not  be  degraded  as  a  result  of  any  orientation. 

3. 2. 5. 8  Nuclear  Hardening.  The  ultimate  roan pack /vehicular  set 
shall  meet  the  nuclear  vulnerability/survivability  criteria  as  set 
forth  in  YEN- 74- 98. 

3.3  Design  and  Construction. 

3.3.1  Materials,  Parts  and  Processes.  Materials,  parts  and  pro¬ 
cesses  used  in  construction  of  the  equipment  shall  conform  where 
practical  to  military  standard  -  consistent  with  achieving  reliable 
performance  under  the  environmental  conditions  3.2 .5. 

3.3.2  Electromagnetic  Radiation.  In  order  to  avoid  undesirable 
interference  effects  on  other  colocated  communications  -  electronics, 
the  equipment  shall  be  designed  to  minimize  electromagnetic  radi¬ 
ation  in  the  frequency  range  of  10  KHz  to  12.6  GHz. 
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3,3  Identification  and  Marking. 

3.3.3. 1  Finish.  Final  paint  shall  be  luster less  green  olive  drab 
as  approved  by  the  procuring  activity. 

3.3.4  Workmanship.  The  equipment  shall  be  manufactured  and 
assembled  in  accordance  with  requirement  9  of  MIL-STD-454. 

3.3.5  Interchangeability.  Parts,  modules,  etc.,  shall  be  inter¬ 
changeable  with  similar  parts,  modules,  etc.,  from  another  set. 

3.3.6  Safety.  The  equipment  shall  be  designed  to  preclude 
safety  hazards  in  accordance  with  the  provisions  of  MIL-STD-454 
Requirement. 

1.  All  corners  and  edges  shall  be  rounded  or  rolled  and  the  battery 
compartment  shall  provide  for  gas  pressure  relief,  as  required. 

3.3.7  Hum^n  Performance /Engineering.  Principles  of  Human 
Engineering  in  accordance  with  MIL- STD- 1472A  shall  be  applied. 

3.4  Documentation  (TBS) 

3.5  Logistics  (TBS) 

3.6  Personnel  and  Training  (TBS) 

3.7  Major  Component  Characteristics.  The  manpack  set  shall 
receive  the  navigation  signals  transmitted  by  the  GPS  transmitters 
and  provide  the  user  with  three-dimensional  position  and  time  in¬ 
formation,  To  accomplish  this,  for  purposes  of  description,  the 
manpack  set  may  bi  considered  as  consisting  of  six  major  functional 
elements  -  an  antenna  assembly,  receiver  assembly,  processor,  com¬ 
puter  programs,  display  unit,  battery  pack,  and  interface  adapters. 

3.7.1  Antenna  Assembly.  The  antenna  assembly  shall  consist  of 
the  antenna,  transmission  lines,  filters,  etc,,  required  to  direct 
the  available  GPS  signals  to  the  rf  input  of  the  receiver  assembly. 
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I 
i 

i  . 

! 

3.7. 1.1  Frequency  Band.  1.558  GHz  to  1.592  GHz. 

3. 7. 1.2  Polarization.  RHCP 

3. 7.1.3  Gain.  Gain  shall  exceed  -1  dbi  over  90%  of  hemisphere 
above  5  degrees  elevation. 

3. 7. 1.4  Antenna  Assembly.  The  manpack  set  shall  contain  an  inter¬ 
nal  antenna  assembly  as  an  integral  part  of  the  package.  Provisions 
shall  be  made  for  attachment  of  a  remote  external  antenna. 

3.7.2  Receiver  Assembly.  The  receiver  assembly  shall  accept  and 
operate  on  the  output  of  the  antenna  assembly.  In  turn,  the  re¬ 
ceiver  outputs  the  required  pseudo- range,  pseudo-range  rate,  and 
system  data  extracted  from  navigation  signals  to  the  data  processor/ 
computer  via  appropriate  interface  circuitry.  Each  receiver  shall 
accept  data  from  the  data  processor /computer  as  required  to  acquire 
and  track  with  the  GPS  signals. 

3.7.3  Data  Processor /Computer  Software.  The  data  processor/ 
computer  software  accepts  and  processes  the  outputs  of  the  receiver 
assembly  and  control/display  units  using  appropriate  hardware, 
software,  and  interfaces  and  outputs  the  proper  quantity,  quality, 
and  form  of  data  required  by  the  remaining  elements  to  support  the 
operations  of  the  set.  As  a  minimum,  the  data  processor /computer 
software  shall  perform  the  following  functions: 

a.  Determine  and  use  the  set  of  GPS  transmitters  from 
those  available  that  will  provide  best  user  position  and/or  navi¬ 
gation  performance  on  a  continual  basis. 

b.  Provide  navigation  signal  acquisition  and  tracking- 
aiding  data  to  the  receiver. 

c.  Accept  control  commands  from  the,  control/display  unit 
and  provide  for  altering  the  processing  of  computer  programs  or 
other  equipment  functions  as  required. 

i 

I 
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d.  Convert  the  system  data  and  pseudo- range /range- rate 
measurements  into  three-dimensional  position  (UTM/UPS,  WGS-72  or 
local  datum)  and  time. 

e.  Provide  for  go /no- no  and  self- test  functions  to  the 
maximum  extent  practical. 

f.  Compute  and  output  to  the  Control/Display  on  command: 

(1)  Present  position 

(2)  Position  of  selected  rendezvous  point 

(3)  Azimuth  to  selected  rendezvous  point  with  refer¬ 
ence  to  true  North 

(4)  Distance  to  selected  rendezvous  point 

(5)  Time-of-day  (day-of-week,  hour-min-sec) 

g.  Provide  as  a  continuous  output  to  the  control/display, 
the  number  of  GPS  transmitters  being  tracked. 

h.  Compute  and  output  on  command  the  CEP/PE  of  position, 
based  on,  for  example,  a  best  estimate  of  range  error  and  GDOP. 

3.7.4  Control/Display  Unit.  As  a  minimum,  the  control/display 
unit  shall  provide  for  manual  input  and  display  output  functions 
including,  but  not  limited  to  the  following: 

a.  Allow  the  operator  to  manually  input  his  approximate 
position,  time,  current  altitude,  etc.,  as  needed  for  processing 
in  accordance  with  3.4.3. 

b.  Allow  the  operator  to  manually  input  up  to  (TBD) 
rendezvous  coordinates. 

c.  Allow  the  operator  to  display  position  and  time  infor¬ 
mation  either  continuously  or  on  command  at  his  option. 

d.  On  command,  display  selected  rendezvous  locations. 

e.  On  command,  display  azimuth  to  selected  rendezvous 
locations,  with  reference  to  true  North. 

f.  On  command,  display  distance  to  selected  rendezvous 

points. 

g.  On  command,  display  the  number  of  GPS  transmitters 
being  tracked. 

h.  On  command,  display  the  CEP/PB  of  position  in  meters. 
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i .  Display.  A  control/display  shall  be  provided  as  an 
input/output  facility  for  the  operator.  As  a  minimum,  it  shall 
permit  day/night  viewing  of  an  alpha-numeric  display  of  worldwide 
horizontal  position  in  either  12-character  UTM/UPS,  or  in  WGS  or 
local  datum  LAT/LONG  coordinates,  plus  4-digit  altitude  in  meters 
relative  to  mean  sea  level  and  system  time  or  local  standard  time 
at  the  option  of  the  user. 

j .  Warning  Lights.  Warning  lights  shall  be  provided  to 
alter  the  operator  to  undesirable  operating  conditions.  These  may 
include,  for  example,  loss  of  carrier  or  code  track  for  each  satel¬ 
lite,  large  GDOP,  poor  data  threshold  and  large  time  bias.  Where 
such  lights  are  provided,  they  shall  be  normally  off,  automatically 
activated,  and  shall  blink  when  activated. 

3.7.5  Battery  Pack.  The  manpack  will  be  powered  by  a  military 
standard  12  volt  or  24  volt  battery  pack  (such  as  the  NICAD  Types 
88585,  88590,  or  Lithium  Types  8A5585,  8A5590  all  of  which  are  now 
under  development).  To  this  end,  the  contractor  shall  provide 
designs  compatible  with  the  use  of  a  battery  assigned  by  the  Power 
Sources  Technical  Area,  ECOM.  Specific  battery  assignment  shall 
be  obtained  by  submitting  AMC  Form  2033- R  through  the  contracting 
officer  to  the  Power  Source  Technical  Area,  ATTN:  AMSEL-TL-P, 

Fort  Monmouth,  N.J.  07703.  The  design  shall  permit  rapid  replace¬ 
ment  of  the  battery  pack,  and/or  recharging  from  an  external  source. 
Reverse  polarity  protection  shall  be  included  in  the  external  power 
interface. 

3.7.6  Interface  Adapter.  The  manpack  design  shall  provide  for 
the  following: 

a.  Communication.  In  order  to  demonstrate  transmission 
of  position  location  and  time  data,  a  digital  output  signal  shall 
be  provided  at  2.4  Kbps  to  modulate  the  AN/PRC-25  or  -77  radio  set. 

b.  Test  Output,  Consistent  with  his  design,  the  contractor 
shall  provide  outputs  useful  for  troubleshooting  and/or  performance 
monitoring/recording  via  a  single  multi- pin  connector.  These  out¬ 
puts  may  include,  but  are  not  limited  to,  such  parameters  as  key 
voltages,  key  waveforms,  pseudo  range /range- rate  information, 
satellite  data,  and  position  and  time  data. 


Page  15  of  15 


c.  Cables /Connectors,  Both  rf  and  power  cables  shall  be 
provided  to  interconnect  with  a  remote  antenna  or  power  source  as 
required. 

d.  Vehicle  Adapter  Unit.  A  vehicle  adapter  unit  shall  be 
provided  to  achieve  the  required  electrical  and  mechanical  inter¬ 
face  associated  with  both  wheeled  and  track  vehicle  use. 

e.  Rack  Adapter  Unit.  An  adapter  unit  shall  be  provided 
to  permit  mounting  the  manpack  in  a  standard  19"  equipment  rack. 

4.0  Quality  Assurance  Provisions  (TBD) 

5.0  Preparation  for  Delivery. 

5.1  Preservation  and  Packaging  for  Domestic  Shipment.  The 
equipment  shall  be  preserved  and  packaged  for  mechanical  and 
physical  protection  in  accordance  with  best  commercial  practice. 

The  equipment  shall  be  adequately  cushioned,  blocked,  and  braced, 
using  suitable  materials  and  containers  as  required  to  afford  max¬ 
imum  protection  from  the  hazards  normally  encourntered  in  storage 
and  transit. 

5.2  Packing  for  Domestic  Shipment.  The  shipment,  packaged  as 
described  in  5.1,  shall  be  packed  in  substantial  commercial  con¬ 
tainers  of  the  type,  size,  and  kind  commonly  used  for  the  purpose 

in  such  a  manner  as  to  afford  maximum  protection  from  normal  hazards 
of  handling  and  transportation  and  to  insure  safe  delivery  and 
acceptance  at  the  designated  point. 


6.0 


Notes. 
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NOTES  ON  BATTERY  CHOICE  FOR  MAil-PACK  EQUIPMENT 
Basic  Types 

All  batteries  can  be  classified  into  one  or  the  other  of  two 
basic  types : 

1.  Primary  -  Non-rechargeable  -  typical  examples: 

Carbon-Zinc,  Mercury,  Silver  Oxide,  Lithium,  Alkaline,  etc. 

*Some  types  of  Alkaline  batteries  can  be  recharged 
20  to  75  times  but  successful  recycling  requires 
care  regarding  overcharging  and  previous  battery 
charge  state. 

2.  Secondary  -  Rechargeable  -  typical  examples: 

Nickel-cadminium,  Lead  Acid,  some  Alkaline  types,  etc. 

Because  of  requirements  related  to  portable  military  field 
equipment,  low  energy  density  primary  batters  can  be  ruled  out 
as  a  feasible  choice.  Mercury  batteries,  while  possessing 
good  energy  density  characteristics,  are  not  suitable  for  opera¬ 
tion  at  temperatures  much  below  0°C  and  therefore  also  are  not 
a  feasible  choice. 

Lithium  batteries  have  an  extremely  high  energy  density 
(highest  of  any  primary  or  secondary  type),  very  good  temperature 
characteristics  and  long  shelf  life.  Because  of  these  desirable 
characteristics,  Lithium  primary  batteries  warrent  further 
consideration  and  are  included  in  the  following  discussion 
along  with  suitable  secondary  battery  types. 
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The  recharageable  Alkaline  battery  as  mentioned  above  is  limited 
to  relatively  few  recharge  cycles  and  further  limited  by  over¬ 
charge/damage  possibilities,  therefore  is  not  considered  suit¬ 
able  for  this  application. 

The  table  summarizes  the  characteristics  of  Lithium,  NICAD  and 
Lead  Acid  (gelled  electrolyte)  batteries,  along  with  other 
factors  considered  important  in  selecting  a  power  source  for 
portable  military  field  equipment. 


Note:  A  sealed  wet  lead  acid  battery  is  now  available. 
(Gates  Rubber) 
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CHARACTERISTIC 


Energy  Density  (W  Hrs/Lb.) 
Nominal  Cell  Voltage  (v) 
Shelf  Half-Life  (at  +25°C) 
Temperature  Range 


LITHIUM 
Primary  Batter 


=150 
2.8  v 


>5  years 


Store 


-60°  to  +74°C 


Operate 

Discharge  Curve  (to  90%  disch) , 
Recharge 
Technique 

Time  (0  to  90%  A  Hr  Rating) 


None 


Life  (#  of  cycles) 
Precautions 

Multiple  Cell  Operation 
Series  (Recom.  Max) 
Parallel  (Recom.  Max) 


No  Restriction 
No  Restriction 


Maintenance  (Bet.  Replacement) 
Mechanical 
Ruggedness 
Package  Style 
Mounting  Position 
State  of  Technology 


Cylinder  cells  &  ? 


New  -  much  dev.  work  in  process 


NIC 

Sealed 


1.2  v 
'=2  months 


-40°  to  +50°C 


-20°  to  +40°C 


SimiL 


Const,  current  (: 


16  hours 


500  -  2000  cycle 
Avoid  low  temp  c 


<8  cells  (9.6  v) 
Not  recommended 


None 


Simil 

Cylinder  cells  ^ 

- Any  ( 

>15  years;  Fast 


Availability  j 1  -  2  sources  -  more  on  way  >5  sources 

Cost  ($/100  W  hr.  Cap.) _ <$20  -  downward  future  trend _ =$100  +  packagi 

(1)  NICAD  cannot  be  recharged  below  +50C;  also  Hz  &  Oz  will  be  released  if  attemp 

(2)  Don't  know  what  recharge  characteristic  is  at  low  temperature  (probably  falls 
temperature  storage  limit  for  discharged  Lead  Acid. 

(3)  Fast  charge  techniques  are  possible  but  are  more  complex  and  subject  to  other 

/  (4)  Hz  &  Oz  is  released  during  charge  (starting  at  =50%  charge)  at  a  rate  of  =8  i 
'  (5)  Multiple  cell  NICAD  should  be  factory-custom  furnished  (cells  must  be  matched 


Multiple  cell  NICAD  should  be  factory-custom  furnished  (cells  must  be  matched 
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LITHIUM 
tr 


NICAD 

(Sealed  Cell) 


~15 
1.2  v 
!=r2  months 


: -40°  to  +50°C 


-20°  to  +40°C  ^ 

-  Similar 


Const,  current  (<A  Hj[/ate)  (3>  (1> 


16  hours 

500  -  2000  cycles 


t 


Avoid  low  temp  charging 


(1) 


<8  cells  (9.6  v)  recommended 
Not  recommended 


(5) 


None 


&  ? 


Similar 
!  Cylinder  cells  ^ 


Any 


(1)?. 


•  more  on  way 
Id  future  trend 


•  work  in  process  >15  years;  Fast  charge  recent 


LEAD/ANTIMONY 

(Gelled  Electy) (Sealed  Ce 11) 

~i5 


2  v 

~16  months 


-60°  to  +60° 
-60  to  +60° 


o)  (4; 

Const.  Voltage  (current  limited) 


16  hours 

200  -  500  cycles 
Gassing^ 


Available  in  2,  4,  6,  &  12  v 
packs  -  no  restriction  on  series 
or  parallel 


Rect.  Packs 

.<*) 


>5  sources 

J£$  100  +  packaging  -  s table 


>2  years  -  some  dev.  work  in 
process 

>2  sources  -  may  be  more  on  way 
=s$30  -  slight  downward  trend 


o  Hz  &  Oz  will  be  released  if  attempted. 

{is  at  low  temperature  (probably  falls  off  fast  below  -20°C.  Also  don't  know  low 
Lead  Acid. 

{are  more  complex  and  subject  to  other  problems.  ) 

p.ng  at  ^50%  charge)  at  a  rate  of  -8  in^  Hz/Hr/100  W  hr.  capacity.  Requires  venting, 
tom  furnished  (cells  must  be  matched  to  minimize  cell  reversal  problems). 
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APPENDIX  C 


A.  1  Recursive  Filter  Equations 

The  solution  to  the  descrete  linear  filtering  problem  with 
recursive  estimation  given  by  Kalman  [Kalman,  1960 ]  assumed 
that  the  evolution  of  a  n-dimensional  linear  multistage 
process,  X^,  as  a  function  of  time,  tj(,  could  be  described  as 


XK  “  *K  XK-1  =  WK 


where  <t> ^  is  the  state  transition  matrix  between  time  t^  and 
tK-l*  ^as  dimensions  of  n  x  n.  represents  the  process 
noise  originating  in  p  <  n  sources  with  zero  mean  and  a  known 
n  x  n  covariance  matrix,  Q^,  with  rank  p.  Kalman  showed  that 
given  a  measurement  set  consisting  of  m  measurements,  ZK, 
taken  at  time  tK>  the  minimum  mean  square  linear  estimate  of 
XK*  XK*  can  obtained  via  a  recursive  calculation. 

The  solution  proceeds  between  time  t^_^  and  as  follows: 

1)  An  estimate  of  the  new  state  at  time  tK 

based  on  the  t  -  is  made. 

K*  1 


CK  “  0K  XK-1 


2)  The  covariance  at  time  t^  is  estimated  as 


p  -  0K  Pr.!  0k  +  Qk 


3)  The  filter  gain  Kp,  is  calculated  as 


KK  “  P!C  HK  ^HK  PK  HK  +  RK? 
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where  Hv  is  assumed  to  be  a  known  m  x  n  measurement  gradient 
matrix  (the  estimated  direction  cosines  for  pseudo  range 
measurements)  and  is  the  covariance  matrix  of  the  measure¬ 
ment  noise  (m  x  m)  which  is  assumed  to  be  uncorrelated  with 

V- 

4)  The  measurement  set,  Zv,  consisting  of  m  measurements 

K  A  A* 

is  used  to  give  an  improved  estimate  of  XK,  XK  , where 


A  » 


LK 


X„  +  Kk  (Zk  -  Hk  Xv) 


K' 


5)  The  covariance  of  Xj.  is  improved  by 


PK=  PK-  KKHKPK 


A  major  problem  with  the  standard  Kalman  is  that  the  matrix 
subtraction  involved  in  the  covariance  calculation  in  Step  5 
can  result  in  a  matrix  which  is  computationally  non-positive 
with  resulting  divergence  of  the  filter.  The  Joseph  form 
of  the  Kalman  filter  modifies  the  covariance  calculation  by 
using  a  less  sensitive  quadradic  form 


K 


(I-KHj^)  PR  (I-KI^)1  +  Kl^  KT 


However  the  Joseph  form  requires  a  considerably  larger  number 
of  computations. 


Square  root  filters  were  proposed  as  a  computationally  more 
efficient  way  to  guarantee  that  is  non-negative  definite 
by  keeping  the  covariance  matrix  in  a  '’square-root"  form  where 
the  square-root,  S^  can  be  obtained  by  Cholesky  decomposition 
of  P^.  The  original  form  proposed  by  Potter  (Battin  1964) 
incorporated  measurements  one  at  a  time  as  scalars.  The 
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covariance  update  is  given  by  a  root -sum- square  (RSS)  procedure 
SK  =  <WK  w*  =  QR)* 
where  WR  -  0R 

Alternatively  the  matrix  can  be  generated  via  a  Householder 
triangularization  procedure  (Kaminski  1971 )  which  is  less 
sensitive  to  computer  word  length  round-off  effects  at  the 
expense  of  more  computation  time.  Kaminski  showed  that  in 
ill-conditioned  problems  square-root  filters  yeild  a  factor 
of  two  more  precision. 

Carlson  ( Carlson, 1973)  proposed  a  modification  to  the  Potter 
form  which  maintains  SK  in  triangular  form..  Thus  in  the 
calculation  of  WK  only  the  elements  in  the  upper  triangular 
part  of  SR  must  be  considered,  reducing  the  computation  by 
one-half.  Carlson  also  gives  a  recursive  calculation  for 
that  does  not  require  matrix  multiplication. 

One  useful  byproduct  of  maintaining  SK  in  upper  triangular 
form  is  that  a  simple  test  for  singularity  is  possible  sine 
will  be  singular  if  one  or  more  of  its  diagonal  elements 
are  zero.  Carlson  shows  that  all  that  is  necessary  to  prevent 
singularity  is  to  require  that  all  of  the  diagonal  elements 
of  Sj(  be  greater  than  the  maximum  round-off  error. 
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A. 2  Orbit  and  Constellation  Selection 

The  purpose  of  this  package  is  to  determine  which  satellites  are 
visible  to  the  user  at  a  particular  time  and  location  and  to 
provide  the  constellation  selection  package  with  direction 
cosines  (and  coarse  ranges).  In  this  computation  a  simple 
circular  orbit  and  a  uniform  gravitational  field  is  assumed. 
Using  orbital  parameters  for  each  of  the  satellites  previously 
stored  in  the  computer  and  the  current  time  the  set  of  satellites 
available  to  the  user  can  be  found  by  a  binary  tree  search. 

The  direction  cosines  to  each  of  the  visible  satellites  is 
passed  to  the  constellation  selection  package. 

The  program  proceeds  in  five  steps: 

1.  Compute  the  inertial  XYZ  coordinates  of  a  candidate 
spacecraft. 

2.  Transform  the  satellite  coordinates  into  earth  fixed 
coordinates. 

3.  Transform  the  earth-fixed  coordinates  into  topocentric 
coordinates  using  the  estimated  user  latitude  and 
longitude. 

4.  Compute  direction  cosine  and  elevation  angle. 

5.  Compute  GDOP. 


* 
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A. 2.1 


Satellite  Motion  in  Orbital  Plane 


Kepler's  equation  is  a  transcendental  equation  in  E(t),  the 
eccentric  anomaly,  given  by 

(t  -  tp)  =  E  (t)  -  e  sin  E  (t) . 

where  P  *  period  of  satellite 

tp  =  time  of  pericenter  in  system  time 
e  =  orbital  eccentricity 
t  =  system  time 

The  first  step  in  the  problem  of  finding  the  satellite  position  in 
user  coordinate  is  to  write  coordinates  of  the  satellite  in  an 
earth-centered  inertial  system  (U,  V,  W)  where  U  is  defined  as 
the  unit  vector  in  the  direction  of  intersection  of  the  orbital 
and  equatorial  plane,  W  is  the  unit  vector  perpendicular  to  the 
orbit  plane  and  V  completes  the  right-hand  system. 

The  position  of  the  satellite  in  the  U,  V,  W  system  is  given  by 
U  -  a(co8  E(t)  cos  u  -  a  (1-E^)^  sin  E(t)  sin 

2  k 

V  -  a(cos  E(t)  cos  u  -  a  (1-E  )’  sin  E(t)  cos  w 
W  *  W(t)  (allows  for  any  out  of  plane  motion), 

where  a  is  the  length  of  the  semi-major  axis  and  u>  is  the  argument 
of  the  pericenter. 

For  the  special  case  of  the  circular  orbit  the  eccentricity  is 
equal  to  zero  and  the  eccentric  anomaly  is  given  by 

E(t)  -  (t-tp),  circular  orbit 

where  tp,  the  time  of  pericenter  passage,  is  interpreted  as  the 
time  of  equatorial  plane  passage  since  in  the  case  of  circular 

C-5 


orbits  t  is  undefined.  Defining  t  in  this  manner  allows  a> 

P  P 

to  be  taken  as  zero. 

However,  for  a  circular  orbit  lying  in  the  equatorial  plane 

(corresponding  to  i  =  0)  the  above  definition  of  tp  must  be 

further  modified  since  the  ascending  node  is  no  longer  defined. 

For  this  case  it  is  convenient  to  take  t  as  the  time  that  the 

P 

satellite  orbit  intersects  the  vernal  equinox,  which  allows 
to  be  set  equal  to  zero.  Alternately,  tp  can  be  any  convenient 
time  and  is  the  angle  between  the  satellite  and  the  vernal 
equinox  at  tp. 

Because  the  simple  Keplerian  model  does  not  account  for  the  non- 
uniform  gravitational  field  of  the  earth  and  the  perturbing 
forces  of  the  moon,  sun  and  other  planets  over  a  period  of  time 
the  error  between  the  actual  and  predicted  satellite  position 
will  gradually  increase  until  acquisition  becomes  difficult 
because  of  doppler  prediction  error  as  well  as  possibly  search¬ 
ing  for  a  satellite  that  is  erronously  assumed  to  be  visible. 

The  orbit  model  parameters  must  be  updated  sufficiently  often 
to  prevent  this  occurrence  either  as  operator  or  satellite  data 
inputs  or  by  having  a  very  large  data  table  stored  in  the  user's 
computer. 

To  provide  the  accuracy  for  navigation,  as  opposed  to  acquisition, 
requires  several  orders  of  magnitude  more  precision  in  the  ability 
to  predict  satellite  position.  This  level  of  precision  is  pro¬ 
vided  by  the  fine  orbit  position  model  whose  parameters  are 
transmitted  to  the  user  as  a  part  of  the  downlink  data.  This 
model  is  used  in  conjunction  with  the  updated  system  time 
estimates  to  refine  the  user  position. 
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2.2.2 


Transformation  of  Satellite  Orbit  Model  Into  User 
Coordinates 


The  coordinates  of  the  satellite  in  the  earth-centered  inertial 
must  be  transformed  into  user  coordinates  for  use  in  the  satellite 
selection  and  acquisition  process  and  in  the  measurement  matrix 
of  the  Kalman  filter.  The  user  coordinate  system  is  a  right- 
hand  system  with  axis  directed  towards  the  local  zenith,  east 
and  north. 

The  necessary  operations  to  transform  the  satellite  coordinates 
in  the  orbital  plane  system  U,  V,  W  into  topocentric  north-east 
zenith  user  coordinates  can  be  written  in  matrix  notation 
(neglecting  flattening)  as 


N™ 

~-sin  0  cos  (3 

-sin  0  sin  0  cos  i 
+cos  0  sin  i 

sin  0  sin  0  sin  i” 
+cos  0  cos  i 

IT 

E 

5* 

-sin  0 

cos  B  cos  i 

-cos  B  sin  i 

V 

z+r-J 

» 

COS  0  COS  0 

cos  0  sin  0  cos  i 
+sin  0  sin  i 

-cos  0  sin  0  sin  i 
+sin  0  cos  i 

W 

where  (U,  V,  W)  are  the  known  coordinates  of  the  satellite, 

<Z>  -  the  geocentric  latitude 

0  ■  true  sideral  time  (e)  +  east  longitude  of  user  (X) 

-  east  longitude  of  the  ascending  node  (ft>) 

i  «  inclination  of  the  orbital  plane 

r  «  distance  from  earth  center  to  user 
o 

N,  E,  Z-  coordinates  of  satellites  in  north,  east,  zenith 

The  elevation  angle  (v)  and  azimuth  (?)  to  the  satellite  in 
topocentric  coordinates  is  obtained  from  the  (N,  E,  Z)  coordinates 
by  y  -  sin  J 

T  -  8l» 

Lcosyj 

Z  E  N 

where  -  »  -  and  -  are  direction  cosines  to  the  satellite 
rs  rs  rs 
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and  r  is  the  distance  to  the  satellite 
s 

(=7n2  +  E2  +  Z2) 

A. 2. 3  Constelation  Selection 

If  more  than  four  satellites  are  visible  to  the  user,  then  the 
user  must  have  a  procedure  for  selecting  the  constellation  which 
produces  the  best  fix.  The  usual  way  of  accomplishing  this  is 
to  calculate  the  Geometric  Dilution  of  Precision  (GDOP)  by  means 
of  the  predicted  direction  cosines  to  each  of  the  batellites. 


Let  1^ 

ml 

nl 

-f 

G  12 

m2 

n2 

-1 

u  - 

X3 

m3 

n3 

-1 

-l4 

m4 

nA 

-1. 

where  1^,  m^  are  the  direction  cosines  from  the  user  to  the  ,:i"  th 
satellite  computed  for  topocentric  coordinates. 

Then  the  GDOP  matrix  corresponding  to  this  particular  constellation 
is  or  11  0-12  o-13  crl4~ 

cr21  (r  22  cr  23  o-24  T  1-1 

E"  <r  31  0-32  <r33  a-  34  “  Gu  Gu 

<r41  0*42  or  43  (r44_ 

i  i  i 

and  Tr  E  |  *  {crll  +  cr2?.  +  cr33  +  cr44}  is  the  GDOP 

Additional  measures  of  position  accuracy  including  circular  error 
probability  (CEP)  in  the  user's  horizontal  poane  can  be  obtained 
from  the  diagonal  terms;  however  the  quantity  of  interest  here 
is  the  GDOP. 

Although  in  principal  the  user's  computer  could  calculate  the  GDOP 
for  each  possible  constellation  of  four  satellites  out  of  the  n 
satellites  in  view,  for  conditions  where  n  is  as  large  as  9  or  10, 
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an  exhaustive  computation  of  all  possible  combinations  would  be 
extremely  lengthy  because  of  the  matrix  inversion  required. 

A  much  shorter  method  which  promises  to  achieve  nearly  equivalent 
results  with  much  less  computation  involves  computing  only  the 
determinant  of  |gu|  and  selecting  the  constellation  which  has 
the  largest  value.  The  reason  this  simpler  technique  can  be 
expected  to  yield  satisfactory  GDOP  discrimination  is  that  large 
GDOP's  result  when  a  particular  coordinate  of  the  topocentric 
system  is  weakly  representated  in  the  system  made  up  of  the 
satellite  position  vectors,  which  leads  to  a  small  eigenvalue 
along  the  weakly  represented  coordinate.  Since  all  of  the 
vectors  which  make  up  the  rows  of  Gu  are  of  the  same  length  a 
small  eigenvalue  will  result  in  a  small  value  for  the  determinant, 
because  the  determinant  of  a  matrix  is  equal  to  the  product  of 
its  eigenvalues.  Although  it  is  fairly  obvious  that  this 
technique  can  be  used  to  eliminate  extremely  poor  constellations, 
it  is  not  clear  that  it  provide  adequate  discrimination  between 
constellations  that  both  provide  a  small  GDOP.  A  short  program 
to  verify  this  procedure  is  being  written. 

Using  one  of  the  above  procedures  the  constellations  are  ranked 
in  order  of  geometrical  performance;  however  several  other 
conditions  must  be  met  before  a  particular  constellation  is 
selected  for  navigation.  One  of  these  is  that  none  of  the 
satellites  in  the  constellation  should  be  below  5°  to  10°  in 
elevation  (unless  it  is  the  only  constellation  or  all  of  the 
other  constellations  have  large  GDOP) .  This  condition  prevents 
excessive  bias  errors  caused  by  ionospheric  and  tropospheric 
delays.  Another  desirable  feature  of  the  selected  constellation 
is  that  none  of  the  satellites  should  be  setting  below  the  low 
elevation  limit  before  the  next  constellation  update  (on  the 
order  of  fifteen  -  thirty  minutes) . 

The  conditions  can  be  incorporated  into  the  basic  GDOP  program 
by  multiplying  the  calculated  GDOP  by  a  set  of  weighting  factors 
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(based  on  experience)  which  will  cause  constellations  with  low 
elevation  angle  satellites  to  hr,ve  a  low  rank  (large  weighted 
GDOP) .  Thus  these  relatively  poor  constellations  will  not  be 
chosen  unless  the  remaining  constellations  have  extremely  poor 
geometry. 

Another  small  point  worth  consideration  is  the  possibility  th§t 
the  selected  constellation  could  become  planar,  causing  Gu  G! 
to  be  singular,  before  the  next  constellation  update.  This 
situation  can  be  predicted  by  examining  the  GDOP  for  the  next 
update  to  see  if  the  determinant  will  change  sign.  If  it  does 
it  will  probably  be  sometime  in  the  interval  where  the  GDOP 
is  very  poor.  Selection  of  this  constellation  should  be  avoided 
whenever  possible. 
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A. 3  Ionospheric  and  Tropospheric  Correction  Techniques 

The  effect  of  the  ionosphere  and  troposphere  is  to  increase  the 
absolute  delay  of  propagating  signals  resulting  in  a  measured 
range  which  is  greater  than  the  actual  range. 

Ionospheric  delay  can  be  corrected  by  a  two  frequency  measure¬ 
ment  or,  with  less  accuracy,  by  predictive  models  of  varying 
degrees  of  sophistication.  Models  avoid  the  additional  hardware 
required  by  the  two  frequency  correction  but  require  good  know¬ 
ledge  of  ionospheric  parameters  for  precise  range  measurements. 

With  a  two  frequency  measurement,  the  range  correction  is 

ar  =  [Pl  -  (f^)2  p2]y[i  -  (f2/fi>2] 

where  p^  and  P2  are  the  measured  ranges  at  frequencies  f^  and 
f2  respectively  (f^>f2>.  This  formula  is  only  valid  at  fre¬ 
quencies  above  1  GHz.  Note  that  all  of  the  necessary  parameters 
are  available  directly  from  receiver  measurements. 

If  two-frequency  measurements  are  not  available  the  ionospheric 
range  correction  can  still  be  determined  using  values  predicted 
by  ionospheric  models. 


One  such  model  |tRW,  1974 |  gives  the  range  correction  R  (in 
feet)  as 


U  •  1.32  (1  +  jf  )  I/f2  <sln2E  1-  §5  +  %)111 

e  e  R_ 

e 


V7here  H  88  ionospheric  scale  height 
Re  ■  radius  of  earth  (6378  km) 

Iy  *  vertical  election  content  (hexams) 
E  88  elevation  angle 
f  -  measurement  frequency  in  GHz 
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Use  of  this  equation  requires  good  knowledge  of  the  ionospheric 
parameters  H  and  I  ,  particularly  at  low  elevation  angles.  If 
these  parameters  are  not  available  to  the  user  it  is  still  possible 
to  make  an  approximate  correction  by  substituting  average  values 
of  H  =  350  km  and  I  =  50  hexams  (the  median  world-wide  daily 
peak)  into  the  above  formula,  where  E  is  known  from  orbital  pre¬ 
dictions.  Under  extreme  conditions  where  I  is  quite  large  or 
changing  rapidly  (as  in  the  polar  cap  region)  the  accuracy  of 
the  correction  using  this  formula  using  measured  ionospheric 
parameters  is  probably  no  better  than  several  meters  and  is  even 
worse  if  average  values  are  used. 

The  two  frequency  measurement  is  capable  of  providing  an  order 
of  magnitude  better  accuracy.  The  primary  restriction  is  that 
the  measurements  should  be  nearly  simultaneous  so  that  the  ray 
paths  traverse  essentially  identical  ionospheric  election 
densities. 

Unlike  the  ionospheric  delay,  tropospheric  delay  is  not  correc¬ 
table  by  the  two  frequency  technique  since  the  group  refractive 
index  is  an  essentially  constant  function  of  frequency  for  radio 
frequencies.  Thus  the  delay  correction  must  be  made  using  models. 
Neglecting  the  tropospheric  delay  would  cause  a  22,4  m  range 
error  for  a  5.7°  elevation  angle  (0.1  radian)  and  11.9  m  error 
for  an  11.4°  elevation  angle  (0.2  radians)  [Thayer,  1967J .  The 
path  length  increase  due  to  ray  bending  can  be  neglected  since 
it  is  less  than  0.5  m  above  5.7°  elevation. 

An  approximate  correction  formula  suitable  for  small  computers 
which  includes  both  dry  and  wet  terms  is  £ Moffett,  1973 J 

ARTR0P0  =  KdP  C8C  [e2  +  ed]  *  +  K*  csc  +  ew 
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where  =  2.278  X  10"^  m/ millibar 

P  =  pressure  at  antenna  height  (millibars) 

E  =  e  lev/at  ion  angle 

V2-5° 

ew"1*5° 

and  is  a  function  of  latitude,  season  and  weather.  Typical 
values  of  are  0.28  m  in  the  tropics  or  midlatitude  summer, 

0.20  m  for  midlatitude  spring  or  fall,  0.12  m  for  midlatitude 
winter,  and  0.05  m  in  the  polar  regions.  A  slightly  simpler 
correction  formula  which  combines  wet  and  dtry  terms  assuming  a 
50%  RH  and  288°K  temperature  is  £  TRW,  1974j 

ARTR0P0  =  7.8  CSC  (E)  exp  |-h/2300C>] 

where  h  is  the  user  altitude  in  feet. 

In  the  above  expressions  the  dry  term  is  an  order  of  magnitude 
larger  than  the  wet  term  for  elevation  angles  greater  than  5°. 

Thus  for  low  accuracy  users  there  is  no  requirement  for  tempera¬ 
ture  or  humidity  inpjts;  however  either  pressure  or  at  least 
height  is  necessary  (which  can  be  converted  to  pressure  assuming 
a  standard  atmosphere).  A  10%  error  in  pressure  will  result  in 
a  2.5  meter  range  error  for  E  =  5°  and  a  1.25  m  error  for  E  e  10°. 


C-13 


mgg 


Jric  I  ass  1 1 1 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whin  Dali  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.I  3-  RECIPIENT’S  CATALOG  NUMBER 


4.  TITLE  (and  Subtltla) 


NAVSTAR  Global  Positioning  System  Special  Studies 
and  Engineering  Program,  Volume  III 


7.  AUTHOROO 

James  J.  Spilker,  Jr. 
Francis  D.  Natali 
Jackson  T.  Witherspoon 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Stanford  Telecommunications  Inc. 
1161  San  Antonio  Road  L 

Mountain  View  CA  94043 


1  I.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

HQ  SAMSO/YEE 

P0  Box  92960,  Worldway  Postal  Center 

Los  Angeles  CA  90009  _ 


4.  MONITORING  AGENCY  NAME  4  AODRESSfif  different  from  Controlling  Office) 


5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final  Report 

March  1974  -  October  1975 


6.  PERFORMING  ORG.  REPORT  NUMBER 

STI  TR  10255,  Vol  III 


8.  CONTRACT  OR  GRANT  NUMBERfa) 

F04701-74-C-0310 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 

NAVSTAR  Global  Positioning 
System  (GPS) 


12.  REPORT  DATE 

October  1975 


13.  NUMBER  OF  PAGES 

148 


15.  SECURITY  CLASS,  (of  thie  report ) 

Unclassified 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

N/A 


16.  DISTRIBUTION  STATEMENT  (Of  thle  Report) 


DISTRIBUTION  STATEMENT  B:  Distribution  limited  to  U.S.  Government  agencies  only; 
Test  and  Evaluation:  October  1975.  Other  requests  for  this  document  must  be 
referred  to  SAMSO(YEE). 


17.  DISTRIBUTION  STATEMENT  (of  the  abetract  entered  In  Block  20,  II  different  from  Report) 


19.  KEY  WORDS  (Continue  on  reverse  aide  If  neceaaary  and  Identify  by  block  number) 


I 


GPS 

manpack  receiver 
direct  acquisition 
synchronization  (GPS) 
table  oscillators 


0.  ABSTRACT  (Continue  on  reverie  aide  If  neceaaary  and  Identify  by  block  number) 


This  volume  presents  the  results  of  the  Manpack  Study  and  Direct  Acquisition 
Study  performed  under  this  contract.  This  report  looks  at  the  feasibility  of 
directly  acquiring  the  GPS  Precision  signal  directly  without  the  aid  of  the 
Clear  signal.  This  report  also  addresses  the  design  criteria  and  specification 
of  a  man-portable,  or  manpack  GPS  receiver.  * 


FORM 
I  JAN  73 


EDITION  OF  1  NOV  85  IS  OBSOLETE 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


